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Multiplicity of compact group representations 
and applications to Kronecker coefficients 

Velleda Baldoni, Michele Vergne 


Abstract 

These notes are an expanded version of a talk given by the second 
author. Our main interest is focused on the challenging problem of 
computing Kronecker coefficients. We decided, at the beginning, to 
take a very general approach to the problem of studying multiplic¬ 
ity functions, and we survey the various aspects of the theory that 
comes into play, giving a detailed bibliography to orient the reader. 
Nonetheless the main general theorems involving multiplicities func¬ 
tions (convexity, quasi-polynomial behavior, Jeffrey-Kirwan residues) 
are stated without proofs. Then, we present in detail our approach 
to the computational problem, giving explicit formulae, and outlining 
an algorithm that calculate many interesting examples, some of which 
appear in the literature also in connection with Hilbert series. 
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1 Introduction 


Let A be a compact connected Lie group and let K be the set of classes of 
irreducible finite dimensional representations of K. Let R be a representation 
of A, which is a direct sum (possibly infinite) of irreducible finite dimensional 
representations of A with finite multiplicities. We write 

V = e^^j^rriK (7r)Vf 

where is the irreducible representation of A parameterized by vr. We 
call the function tt —)■ mi^(7r) on A the multiplicity function of the repre¬ 
sentation V. The study of the multiplicity function m/f(7r) is important in 
representation theory, invariant theory, quantum information theory. 
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When the representation space V is constructed by ” geometric quanti¬ 
zation” of a Hamiltonian manifold M, the moment map <h on M gives us 
a geometric interpretation of the multiplicity function of the representation 
V. This is the famous [Q, i?] = 0 theorem, obtained by Meinrenken-Sjamaar 
[TT] , This has consequences on the qualitative properties of the function 
TT —)■ m;^(7r) that we will recall here in two examples, which are the paradigms 
of geometric quantization: 

• The quantization of a symplectic vector space under a linear action of 
K. 

• The quantization of T*K, the cotangent bundle of K. 

In these two basic examples, we give a direct construction of the corre¬ 
sponding representation space V. We will not justify that the representation 
space V is “the quantized space” of M, but we will recall results on the 
multiplicity function of the representation of K in V in terms of the moment 
map on M which conhrm the fact that V is the “the quantized space” of the 
Hamiltonian space M. 

Example 1. Quantization of a symplectic vector space 

Consider M a 2n dimensional symplectic vector space (M = with 
a linear symplectic action of the compact group K. We choose a Hermitian 
structure on M (so we identify M to C") such that K acts unitarily on M. 
We denote by Sym^{M) the complex space of fc-symmetric tensors on the 
complex space M. Consider 

the space of polynomial functions on M*. 

The space V may be completed as a Hilbert space, by choosing a Gaussian 
inner product, and the completion of V is the Fock space. In the philosophy 
of quantization, this Fock space is the quantized space associated to the dual 
representation of K in the symplectic space M*. 

If K contains the homotheties, we then have 

V = 

where mi^(7r) is finite. When K is abelian, mx(vr) is a partition function. 

Remark that the knowledge of the function tt —)■ mx(vr), for vr trivial on 
the semi-simple part of K, allows to compute the Hilbert series of the ring 
of invariant polynomials on M* under the semi-simple group [K,K]. 
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Example 2. Quantization of the cotangent space T*G 

Consider a compact Lie group G. Let us define V = R{G) to be the 
subspace of C°°(G) generated by the coefficients {gui,U 2 ) of finite dimen¬ 
sional representations of G. The space L‘^{G) is the Hilbert completion of V. 
In the philosophy of quantization, the space L^{G) “is” the quantized space 
associated to the symplectic space T*G. 

By Peter-Weyl theorem, under the action of G x G by left and right 
translations: 

Let K be a subgroup of G, and consider the subgroup G x JL of G x G. 
Under the action of G x iC, 

^ ® (U)- 

where p) is the multiplicity of the representation jj, in the restriction 

of A to K. This is the so called branching coefficient. 

In both of these examples, we will recall some of the qualitative proper¬ 
ties of the multiplicity function, in particular its piecewise quasi-polynomial 
behavior on closed cones. These examples being particular cases of geometric 
quantization of Hamiltonian manifolds, the moment cone, its decomposition 
in cones of quasi-polynomiality, and the determination of its faces, gives us 
already information on the behavior of the multiplicity function. 

We will then recall our method to compute the multiplicity function 
m;^( 7 r), based on computations of partition functions via multi-dimensional 
residues. We apply it to the challenging example of computation of Kronecker 
coefficients. In particular, we obtain an algorithm which enables us to com¬ 
pute the dilated Kronecker coefficients associated to Young diagrams with 
n rows (but “ any number” of columns and any shape) in polynomial time 
when n is fixed. We implemented our algorithm as a simple Maple program. 
In particular, it computes the dilated Kronecker coefficient g{ka,k(3,k'y) in 
reasonable time (reasonable=less than 20 minutes) for 3 rows. We obtain a 
quasi-polynomial of degree 11 for general a,/ 3,7 (with coefficients periodic 
functions of period at most 12) see Example 1831 When a, (3, 7 are special, the 
degree might be much smaller. For example, the dilated Kronecker coefficient 
m(fc) = g{k[l, 1,1], k[l, 1,1], k[l, 1,1]) corresponds to the Hilbert series of the 
ring of invariants of SL{3) x SL{3) x SL{3) in 0 0 C^. An efficient 
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way to represent periodic functions is by using step-polynomials ([SI])- In 
this representation, m{k) is given by the following quasi-polynomial 
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Here for s e R, the function {s} = s — [sj G [0,1) where [sj denotes the 
largest integer smaller or equal to s. It is easy to check that (fortunately) 
this result of our algorithm agrees with Kac’s determination [26] of the ring 
of invariants [Sym{C^ <8) (8) C^)]‘S'-f'(3)xSL(3)x5L(3)^ which is freely generated 

with generators in degree 2,3,4 (see [IQ]). Indeed 


k 


1 

(l-t2)(l-t3)(l-t4)- 


Another example is 

m{k) = ^(fc[3,3,3,3],fc[4,4,4],fc[4,4,4]) = dim[^|/m^2fc(c^^c3(8C3)]®^(4)xSL(3)x5L(3)_ 


Our algorithm gives the Hilbert series 


k 


1+t^ 

{1 - {1 - 1 ^) {1 - 1 ) {1 - t^y 


We will resume with the study of Kronecker coefficients in Section 16.21 
giving all the details. For example, we can compute easily the Hilbert series 
for 3,4, 5 qubits. The result for the case of 5 qbits appears also in [38] with 
a discrepancy due to a small misprint fsee 16.2p . 

Our article, focused towards computational problems, is strongly inspired 
by the article of Christandl-Doran-Walter [13] and we thank them for intro¬ 
ducing us to the subject of quantum computing. Our computational method 
differ at some points from their method, and we will discuss the differences. 
It is based on Jeffrey-Kirwan residues. This allows us to produce a symbolic 
function valid inside the cone of quasi-polynomiality, in particular along a 
line. 


We believe that this survey might be of interest. The eventual reader that 
just want to read a discussion of general results in Hamiltonian geometry and 
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some applications of these results to combinatorics, without being addicted 
to computing, might skip the technical parts, namely. Subsection 14.91 inside 
Section H] as well as Subsection 15.21 and Subsection 15.31 inside Section [5l and 
might go directly to Section [ 6 ] skiping again the technical subsection 16.11 

2 Notations 

Let iL be a compact connected Lie group and Tk a maximal torus of K. We 
denote by i and h the corresponding Lie algebras. Denote by Wf the Weyl 
group, and w —)■ e{w) = deti^ w its sign representation. 

The weight lattice of T is a lattice in it^. If A G A^ C ztf , it determines 
a one dimensional representation of by f —)■ with t = exp X, X G t. 

As A takes imaginary values on h, is of modulus 1. 

We denote by Lj^ G iit the dual lattice of Li^: if A G Ak, 7 G L^, then 
(A, 7 ) is an integer. 

Let Af C it* be the root system for t with respect to h. If a G A^, its 
coroot Ha is in ih, and (a, Ha) = 2. Dehne K to be the set of irreducible 
hnite dimensional classes of complex representations of K. Fix Aj^, a positive 
system for A{, and let 

= >0,aG A+} 

be the corresponding positive Weyl chamber. Let pt = \ X]agA+ 

We denote by Ak,>o the ’’cone” of dominant weights, that is the set 
Ak n > 0 - An element A G Ak,>o is called a dominant weight. 

When the group K is understood, we denote Ak simply by A, Tk by T, 
h by t, etc. 

The following example is the only example which will be needed when 
discussing Kronecker coefficients. 

Example 3. 

We consider the case K = U{n) and T G K the torus consisting of the 
diagonal matrices. Then the Lie algebra t consists of the nxn anti Hermitian 
matrices and it is the space of Hermitian matrices. If we identify t and t* 
via the linear form Tr{AB), then t = t* is the set of diagonal anti hermitian 
matrices. Thus the positive Weyl chamber is = {■C = [■Ci;'^ 2 , • • •, ■Cn]} 
with G M and > ^2 ^ ^ where ^ represents the hermitian matrix 
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with diagonal entries ^j. Denote by |^| = the trace of the matrix 

corresponding to 

An element X G H* is conjugated to a unique element ^ ^ 2 , ■ ■ ■, ■Cn] 

in it>Q, where the list .^ 1 , ^ 2 , ■ ■ ■ ,^n is the list of eigenvalues (with their mul¬ 
tiplicities) of X ordered decreasingly. We denote ^ by Spec(X). 

The “cone” of dominant weights is ^u(n),>o = {A = [Ai, A 2 ,..., A„]} with 
Xj G Z and A^ ^ A 2 ^ ^ X^. 

If A G A[/(„)^>o is such that Xn > 0, it indexes a finite dimensional irre¬ 
ducible polynomial representation of GL{n, C). The corresponding subset of 
Af/(„)^>o will be denoted by PA; 7 („)^>o. If A G PA^(„) > 0 , we also identify A 
to a Young diagram with n rows. The content of the corresponding diagram 
is the number of its boxes, that is A; = |A|. The dominant weight [k,k,, k] 
correspond to a rectangular Young diagram with k columns and n rows, and 
index the one-dimensional representation det( 5 ')^ of U{n). 

Assume now N > n, then there is a natural injection from PAu(^n),>o to 
PA[/(Ar),>o obtained just by adding more zeros on the right of the sequence 
A. We denote by A the new sequence so obtained. □ 

We parameterize K by the set of elements (tta, where tta denotes the 
irreducible hnite dimensional representation of highest weight A G Ak,>o and 
is the finite dimensional space on which tta acts. If there is no ambiguity 
on the group K, we may write simply Va- If T is a torus, we also write for 
the one dimensional representation of T associated to a weight of T. The 
contragredient representation of K is indexed by the dominant weight 
A* = —tco(A), where wq is the longest element in >V{. 

Let TT be a representation of K acting on a complex vector space V. As¬ 
sume Y is a direct sum (possibly inhnite) of irreducible finite dimensional 
complex representations tta of K, each of them occurring with finite multi¬ 
plicities, and let 

Assume that the restriction of the representation of K to its maximal torus 
T is also with hnite multiplicities mT{fi). We have 

rriKiX) = ^ e{w)mT{X +Pi-w{pi)). ( 1 ) 

wGWf 

This is a consequence of the denominator formula: 

JJ (1 - e“) = ^ (2) 

oeA+ w€Wt 
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\i X E we denote by K{X) the subgroup of K stabilizing X. If K is 
compact connected, then K{X) is a compact connected subgroup of K. 

If is a Hermitian vector space, dehne Tipure = {v E Ti, {v,v) = 1}, the 
set of elements of H of norm 1 . Such an element is called a pure state. A 
density matrix is a semi-positive dehnite Hermitian matrix with trace 1. A 
pure state in the space is called a qubit, and a pure state in 
called a N qubit. 

For us, a cone C in a real vector space E is a. closed subset of E contain¬ 
ing 0, invariant by positive homotheties. We consider here only polyhedral 
cones, that is ^ G C* if and and only if ^ satishes a certain number of linear 
inequations {Xa,0 > 0, with Xa E E*. When we say cone, it will always 
mean a polyhedral cone. Most of the time, E will be equipped with a lattice 
L. Consider the dual lattice L* in E*. We say that C* is a rational cone if the 
Xa are in L*. In other words, ii E = M"", and L = Z”, a rational cone is a 
cone dehned by inequations with integral coefficients. 

We will say that a cone C is solid if C has non empty interior. If F is a 
hnite set of vectors in E, we denote by Cone{E) the cone generated by E. 

An affine cone is a translate s -I- C of a cone. An open cone will mean 
the interior of a polyhedral cone, that is a set determined by strict linear 
inequations {Xa,0 >0. If C is a polyhedral cone in a vector space, we say 
that C = UCa is a cone decomposition of C if the Ca are closed polyhedral 
cones of dimension equal to the dimension of C, and if, when a ^ b, the 
intersection Ca n Cf, is contained in the boundary of the cone Ca and C^. 

A face F of a cone C generates a linear space that we call lin{E). A facet 
of F is a face of co dimension 1 in F. 

A polytope is a closed compact convex subset of F determined by linear 
affine inequations {Xa,0 ^ C, with Xa E E*, Ca E M. If Xa E L*, and 
Ca E Q, the polytope is rational. 

3 Basic examples for multiplicities 

Example 4. The Littlewood-Richardson coefficients 

We consider two irreducible representations (tta, 14) and (tt^, I 4 ) of K. 
Their tensor product V = V\ ® Vp is a K x K representation with action 
dehned by: (tta 07 r^)(/ci, ^ 2 ) = T^xikf) 07 r^(/c 2 )- The group K acts diagonally 


on V via 7r(k) = TT\{k) ® T^u{k), so we may consider the representation vr = 
tta ® restricted to K. We can write the classical formula 

(Va ® — ®!yeA>oCA,/_i^ 

where is the multiplicity of tIu in (tta ® 7r^)|A'- 

The numbers are called Littlewood- Richardson coefficients. 

The function k —)■ ^ {0; 1)2,...} is called the dilated Littlewood- 

Richardson coefficient. It follows (see Theorem l63|) from the [Q,R] = 0 the¬ 
orem that the function k —)■ is given by a quasi-polynomial formula for 

k>l (polynomial in the case of t/(n), see Proposition l30|) h 

We recall that Cochet, [14] , [15], has given an algorithm to compute the 
dilated Littlewood-Richardson coefficient for all classical root systems. This 
algorithm is available on [53] □ 

Example 5. The Kronecker coefficients 

Let N = nin 2 ■ ■ - ng where ni, 77 , 2 ,..., are positive integers and write 
® ■ ■ ■ ® C"'^ Consider the action of the group K = U{ni) x 
• • • X Uijis) on the complex vector space C^, where (/ci, k 2 ,... ,ks) G iC acts 
hj ki®k 2 ®- ■ -^kg on ®0• • • 0. Consider the space V = , 

where = Sym^{C^) is the space of symmetric tensors of degree k. Write 

V = @j(Ai, A,,..., A.)©'"‘’ ® • ® V®->. 

Here Aj are polynomial representations of U{nj) and are indexed by Young 
diagrams with nj rows. 

Considering the action of the center, we see that all diagrams Xj occur¬ 
ring in have content k, so that they also index irreducible representa¬ 
tions 7i\. of the symmetric group By Schur duality, ^^(Ai, A 2 ,..., A^) 
is the multiplicity of the trivial representation of in ttai 0 ••• 0 t^Xs- 
The numbers 5 '(Ai, A 2 ,..., A^) are called Kronecker coefficients. The function 
k —)■ g{k\i, kX 2 ,..., kXg) of k E {0,1,2,...} is called the dilated Kronecker 
coefficient. It follows again from the [Q, R] = 0 theorem that the function 
k —)■ g{kXi, kX 2 ,..., kXg) is given by a quasi-polynomial formula for k > 1, 
(Proposition [3T|) . 

Denote by M = 77.2773 ■ ■ - Ug. In computing Kronecker coefficients, we may 
assume 77,1 < M, and that 77.1 is the maximum of the n^. Indeed, if ni > M, 
the multiplicities g{Xi, X 2 ,..., A^) stabilize in the sense that g{Xi, X 2 ,..., A^) 
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is non zero only if Ai is obtained from an element in PAu(^m),> 0 j by adding 
more zeroes on its right, and multiplicities coincide. Moreover, we have 

g{\i, A 2 ,..., As) = g{Xi, A 2 ,..., As). 

Thus it is sufficient to study Kronecker coefficients in the case where ni < 
M = n 2 ri 3 ■ ■ - Us and where Ui is the number of rows of the tableau corre¬ 
sponding to Aj. 

To describe an algorithm to compute the dilated Kronecker coefficient is 
the main computational objective of this article. We will state some of the 
previous results and our results in Section El Our algorithm uses a branching 
rule from U{n 2 n 3 ■ ■ -Us) to U{n 2 ) x Uln^) x • • • x U{ns) in order to reduce 
(slightly) the size of the problem. Our Maple program is available on [53] □ 

4 Linear representation of in a Hermitian 
vector space TL 

Let T-Lhe a. finite dimensional Hermitian vector space provided with a repre¬ 
sentation of K by unitary transformations. 

Assume that K contains the subgroup of homotheties Consider 

V = Sym{'H), the space of symmetric tensors, so we have 

Sym{n) = ©AeKm5^(A)K/ 

where m^{\) is finite. 

It is clear that if mj^{a) and m^(/3) are non zero, then m^{a + (3) is 
non zero. Indeed the product of two non zero vectors /, g in Sym{'H) is non 
zero, and if f,g are highest weight vectors of weights a, [3, the product is 
a highest weight vector of weight a + f3. So the support of the multiplicity 
function m^{a) is a “discrete cone” (that is a semi-group). We will relate 
this discrete cone to the moment map and to the Kirwan cone. 

4.1 Examples of decomposition of Sym{'H) 

Example 6. 

Let l-L = £A and let K = acting by the homothety 

{zi, Zn) P^Z2, . . . , P^Zn). 
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Then 


V = Sym{H) 


sOO 

4=0 


fn — 1 + k\ 

V ri-l ) 


^ike 


since dim{Syrn^{'H)) = 

So the mnltiplicity fnnction k —?■ dim{Sym^{'H)) is a polynomial function of 

k. 


Example 7. The Knapsack 

Again, let "H = C” and let K = acting by 

[Zi, Z2, ■ ■ ■ , Zn) ^ [e Zi,e Z2,---,e Zn), 

where now the Ai are any positive integers. 

Then Sym{'H) = '^i.m{k)e^^^ where m{k) is the number of solutions in 
non negative integers x, of the knapsack equation 


AiXi + A 2 X 2 + • • • + AnXn — k. 

The computation of the function m[k) is an ’’intractable” problem, as il¬ 
lustrated in the lecture https://www.youtube.com/watch?v=2IbJf4oX0xU&feature=youtu.be 
by P. Van Hentenryck. See however [3] for results on its highest coefficients. 

Example 8. Cauchy formula 

Let N,n be positive integers, and assume N > n. Let A = [Ai,..., A„] 
be a sequence of weakly decreasing non negative integers with A„ > 0. We 
consider A as an element of Au^n),>o, that is a dominant polynomial weight 
for U{n). To A G Au{n),>o is associated an irreducible representation of U{n) 
that we have denoted by Recall that if iV > n, then there is a natural 

injection A ^ A from PAtj(^n),>o to PA[ 7 (jv),>o obtained just by adding more 
zeros on the right of the sequence A. The decomposition of Sym{CA (g) C'^) 
with respect to U{n) x U{N) (see [36], page 63) is given by Cauchy formula: 

SymiC ® C~) = 0 (3) 


Example 9. Clebsch-Cordan coefficients 
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Consider the representation 11 of Jt" = U{d) x U{d) x U{d) on "H = 
0 [(d) © 0 [(d) given by 

U{g,h, k){A,B) = {gAk-\hBk-^), 

where 0 [(d) is the Hilbert space of complex d x d matrices eqnipped with 
the trace inner product {A,B) := TtAB*. This action factorizes through 
the center Z = oi U{d) embedded in 17(d) x 17(d) x 17(d) as a diagonal 
subgroup {zld, zld, zld). We can write 

Sym{H) = © H;. 

Here A, /i, z/ vary in >o and is the multiplicity of the representa¬ 
tion W in Va © For to be non zero, we need |A| + |p| = \v\, as seen 

by considering the action of the center. 

4.2 The moment cone 

In the following, the compact group K will be fixed, and we denote simply 
by T its maximal torus, t its Lie algebra, etc., as we stated in the Section [2] 
Consider the moment map "H —)■ it* given by 

$i^(n)(X) = (Xn,n). 

Here X E t, and we have denoted by n —)■ Xv the inhnitesimal action of X G ^ 
in 1/ by a anti-hermitian transformation, so {Xv,v) is purely imaginary. 

Remark 10. 

We consider "H as a symplectic manifold, with symplectic form ^(dn, dv) 
(if l-L = with coordinates Zk, this is the form ^J2k^^kdzk), then 
is the moment map for the action of 77 in "H in the sense of Hamiltonian 
geometry. 

□ 

We consider H>o as a subset of it*. 

Definition 11. Define 

.CK{n) = d>K{n)ntt*^, 
and 

• Xk^H) = ^K^kipure) C it>o. 
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Assume that K contains the homotheties and let J G it so that the 

inhnitesimal action of J is the identity on l-i. Then a pure state v is such 

that ^k{v){J) = 1. 

Recall the following theorem, which is a particular case of Kirwan theorem 
[30] (a proof of this theorem, following closely Mumford argument, jH], can 
be found in [S]). 

Theorem 12. • The set CxiTL) = fl it>Q is a rational polyhedral 

cone. We call CxiTL) the Kirwan cone. 

• The set Ak{TL) is a rational polytope. We call Ak{'H) the Kirwan 
polytope. 

Thus there exists a hnite number of elements Xa G T such that 

We say that the inequations (Xa,^) > 0 are the inequations of the cone 

CxiTL). We may normalize Xa to be a primitive element in T, the dual 

lattice to A. The set AxiTL) is the intersection of Cxi'H) with the affine 
hyperplane (J, 0 = 1 


CK{n) = R>oAK{n) 

is the cone over the Kirwan polytope. 

It is in general quite difficult to determine the explicit inequations of the 
cone Ck{TL). An algorithm to describe the inequations of this cone, based 
on Ressayre’s notion of dominant pairs, [IB], is given in Vergne-Walter, [B2] . 
Let us hrst give the emblematic example of the Horn cone. 

Example 13. Horn problem 

Consider 3 Hermitian matrices X,Y,Z. Let Spec(W), Spec(y), Spec(Z) 
denote the list of eigenvalues of the Hermitian matrices X, Y, Z. The Horn in¬ 
equations describe the range of the triple (Spec(W), Spec(y), Spec(Z)), when 
the 3 matrices X, Y, Z are constrained by the relation X + Y + Z = 0. This 
problem is related to the moment map for Example [9] as follows. 

Consider the representation oi K = U{d) x U{d) x U{d) on TL = gl(d) © 
The moment map is 

$x(A, B) = {AA\ BB*, -AM - B*B). 
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Since any non-negative Hermitian matrix can be written in the form A*A 
and since the spectra of AA* and ^4*^4 are equal, the moment cone is equal 
to 

CK(n) := {(Spec(X), Spec(y), Spec(Z)) : X,Y >0,Z <0, X + Y+Z = 0}, 

As proved in [31] , [35] (see also 0 ). the equations of Ck^H) are given by the 
inductive system of inequalities conjectured by Horn [21]. These inequations 
are of the following form. Let /, J, K be subsets of [1, 2,..., d], all three 
of them of cardinal r < d. Let Ei be the diagonal Hermitian matrix with r 
eigenvalues 1 in positions /, and others being 0. Then the triple [Ej, Ej, Ek) 
gives rise to the inequation 

Tr{EjDi) + Tr{EjD 2 ) + Tr{EKDs) < 0 (4) 

on triples {Di, D 2 , D^) of Hermitian diagonal matrices. 

Horn defined inductively, for every r < d, a set H{r, d) of triples (/, J, K) 
of subsets of cardinal r of [1,2, ...,d]. Then Ck{H) is described by the 
inequations dl]) above, for all {I,J,K) G H{r,d) and all r < d, and the 
equation Tr{Di) + Tr{D 2 ) + Tr{Ds) = 0. 

We now consider examples related to Kronecker coefficients. 

Example 14. (Bipartite case) . 

We return to Example (8] Consider H = EA <Z) C'^ with the action of 
U{n) X U{N). Using the Hermitian inner product, we identify A G "H to a 
matrix A : C” —?■ C^. Then the moment map is given by 

$^(A) = [AA*, A*A] 

with value Hermitian matrices of size n and size N respectively. The matrices 
AA* and A*A have the same non zero eigenvalues. Assume that N > n. 
Define 

PKin),>0 = {e = [6, • • • , en]; 6 > 6 > • • • > en > 0}, 

and consider as a subset of the positive Weyl chamber A>o for 

U{n). There is a natural injection from to obtained just 

by adding more zeros on the right of the sequence We denote by the 
new sequence so obtained. Then we see that the Kirwan cone CxiB,) is the 
’’diagonal” (^,|) with ^ G PAl{n),>o- 
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□ 

In the above example, we see that the cone Ck{'H) may have empty 
interior in *t>o- The following general result holds. 

Lemma 15. The cone Ck{H) is a solid cone if and only if there exists v &TL 
so that the stabilizer of v is a finite group. 

Return to the Kronecker case. 

Consider Ti = C"! ( 8 )C"' 2 (g).. .(g)C"'s with action oi K = U (ni) x - ■ - xU (n^) 

on H. 

We may assume ni > n 2 > ■ ■ ■ > Ug, and let M = 77 , 277.3 •••ris. The 
corresponding Kirwan cone Ck{TL) is a subset of 

-P*tu(ni),>0 ® -^*^u(n2),>0 ® ' ' ' ® >q. 

Then if 77i > M, this cone stabilizes in the sense that CxiTL) coincides 
with the cone associated with the sequence (M, 77 . 2 , 773 ,..., Ug) embedded by 
sending to PHh^ni) >0 ^ ^ Furthermore, considering 

the action of the center, we see that Ck{TL) is contained in the subspace 
E = {( 7 / 1 ,..., 7/s)} of itj dehned by the s — 1 linear equalities Tr(7/i) = 
Tr(7/2) = ■ • • = Tr(7/s). 

For later use, we recall the following result (see [S2])- 

Proposition 16. Let TL = <8) • • -(gjC""* with action of K = U (77,1) x 

• • • X Uijig) on Ti. Assume that all Ua are greater or egual to 2, and that s 
is greater or egual to 3. Assume also that ni = max{ni) and 77.1 < 77.2 • • -Ug. 
Then Ck{TL) is solid in the subspace E. 

A few examples of (tti, 77 . 2 ,..., tt.^) where the inequations of the cone 
Ck(TL) are known are: 

(2, 2,..., 2) (with any number N of 2: the cone of iV-qubits,), (Higuchi- 
Sudbery-Sultz p3]h 
(3, 3, 3) (Franz [20]). 

(4.2.2) (Briand [TO], Bravyi [5]), 

(6.3.2) , (9,3,3) (Klyachko [3T]h 

(4,4,4) and (12,3,2,2) (Vergne-Walter, [52]) (this last case being incor¬ 
rect in Klyachko). 

The general case seems for the moment out of reach. 

Example 17. 
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Let us write explicitly Higuchi-Sudbery-Szulc description of the cone of 
A^-qubits (see also [H]). Consider Ai = [Aj, A^],..., Aat = [\^, A^] a sequence 
of N elements of >o ( assume that Tr(Ai) = 

Tr(A 2 ) = • • ■ = Tr(AAr). Then (Ai,..., X^) G C'{/( 2 )x---xc/( 2 )(C^ 0 • • • (8) C^) if 
and only if 


for any j = 1, 2,..., A^. 
□ 


Example 18. Quantum marginals 

Let A,B,C be integers. If n G 8) C® 8 is a pure state, its 
quantum marginals are dehned as follows. Consider v as an operator v : 

—)■ C'® 8 and dehne pa{v) = v*v. Dehne pb{v), pc{v) similarly. Then 
Pa^v), pb{v), pc{v) are densities matrices of size A,B,C respectively. They 
are called the quantum marginals of v. 

Consider "H = 8 C'® 8 C'", with action oi K = U (A) x U(B) x U (C). 

The moment map '■ 8 C-® 8 fu(A)* © ivi{B)* © iu{C)* is 

= {Pa{v),Pb{v),pc{v)). 

Given 3 densities matrices Pa,Pb,Pc, the quantum marginal problem is 
to determine if there exists a pure state v G 8 C'® 8 C*" with quan¬ 
tum marginals pA,pB,Pc- We see that this pure state exists if and only if 
(ypA, Pb, Pc) is in the image of the moment map ^k- It is thus necessary 
and sufficient that the triple of spectrums (Spec(p^), Spec(pB), Spec(pc')) of 
(Pa, Pb, Pc) satisfy a certain number of linear inequalities. Unfortunately, 
there is not a good understanding of what are these inequalities, except in 
the few low dimensional cases quoted previously. 


4.3 Duistermaat-Heckman measure on the moment cone 


For simplicity, we assume that the moment cone CxiB) intersects the interior 
of the Weyl chamber. Of course, this is the case when CkQ-L) is solid. We 
also assume (we can always reduce easily to this case) that the kernel of the 
homomorphism K -A- U{'H) is trivial. 
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Definition 19. Define 


r = dimc('H) — |A^| — dimCi<-('H). 
Then, if the cone Ck{H) is solid, 

r = dinic('H) — |A^| — dimt. 


Consider the Lebesgue measure dv on 1-i (considered as a real vector 
space). If ^ there is a natural measure on the coadjoint orbit 
of S, determined by the symplectic structure of (see H). Consider the 
Duistermaat-Heckman measure supported on the polyhedral cone 

Cxi'H), determined by 


f{^K{v))dv 


'V 



DH«(0. 


Here / is a compactly supported continuous function on it*. In short, we 
divide the push forward {^K)*{dv) of the Lebesgue measure dv on "H by the 
Kirillov measure of the orbits in the image. This measure is strictly positive 
on the relative interior of Ck{'H). 

Consider the Lebesgue measure d^ on the vector space lin{Cxifid)) spanned 
by Ck{'H). Here we normalize d^ so that it gives mass 1 to a fundamental 
domain of A n lin{Ck{'H)). 

The following theorem follows from Duistermaat-Heckman [T9] . 


Theorem 20. There exists a cone decomposition CKiTL) = UaeFCa? m closed 
polyhedral cones Ca, and for each a, there exists a homogeneous polynomial 
function ^ of degree r on lin{CK{H)) such that 




%ffe Ca. 

We denote by dh^ the function on CxifH) = UCq such that dh^ = on 
Ca. So we obtain a piecewise polynomial function dh^ with support CxiTL) 
and continuous on Cxi'H). We call the Ca chambers of polynomiality (for 
the Duistermaat-Heckman measure). It is already quite difficult to determine 
the cone Ci^('H), so even more so to determine the chambers of polynomiality 

The function dh^{f) is related to the volume of the reduced fiber of TL. 
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Definition 21. Let ^ E it*, we define the reduced fiber ofT-L at ^ by Hredi = 

ixim/K- 

In other words, the reduced space at ^ consists in all orbits Kv of K in 
LL projecting on the orbit under the moment map. Via the identification 
of 'Hred,i with the GIT (geometric invariant theory) quotient of semi-stable 
orbits under Kc, the reduced space can be provided with a structure of 
projective variety (see [ 12 ] )• 

If all orbits of K in have the same dimension, the reduced space 

is an orbifold. We say that ^ is a quasi-regular value. In particular, if the 
map k ^ kv \s injective for all n G "H projecting on we say that the action 
of K on is free. In this case ^ is a regular value of <hi^, and Tired,^ 

is a smooth manifold. Furthermore, Tired,^ inherits a symplectic structure 
fig from the symplectic form fl on Ti: the restriction of fl to is 

the pull back of a form fig on Tired,i- For ^ in the relative interior of the 
Kirwan cone, and a quasi-regular value, then the reduced space Tired,^ is 
of dimension r. Duistermaat-Heckman theorem implies that dfifi^ifi) is the 
volume of the symplectic space Tired,i- So when we restrict the function 
to the line t^ (^ in the interior of Gx('H) and quasi-regular), the function 
t -E- dhfi^it^) is homogeneous of degree equal to r. This is also true without 
the quasi-regularity assumption. 

Lemma 22. //^ is in the relative interior of the cone CxiTi), then dhffi-{tf) = 
t^dh'^if). 

So the homogeneous degree on a line stays constant for interior lines. But 
typically, the function dh^ vanishes on the boundary, except if it is constant 
on Gi^('H). 

Example 23. (The Bloch sphere) 

Let us consider a very simple example. 

Let T be a two dimensional torus, with Lie algebra t = MJi © MJ 2 . Let 
t* with dual basis J^. We consider the diagonal action of T on "H = C^, 
by exp(diJi + 6*2 J 2 )('^i, 2 : 2 ) = Zi, Z 2 )■ Thus the weights 0 i ,02 of the 
action of T on are = iJ^, 02 = iJ^. 

The space of pure states divided by the action of (e*®,e*®) is the Bloch 
sphere. The Kirwan polytope is the segment [ 0 i, 02 ] in © 1^02 = it*- 
The Kirwan cone is the cone M>o0i © R>o 02 - The reduced fibers are points, 
and the Duistermaat Heckman measure is the characteristic function on this 
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cone. This is just Archimede result for which the area on the sphere between 
two latitudes depends only of the difference of their heights on the z axes as 
one can see looking at the picture Figured! 


2 





r 


—r 


Figure 1: The Kirwan polytope 


4.4 Moment map and multiplicities 

We continue (for simplicity) to assume that Ck{'H) intersects the interior of 
the Weyl chamber and that the kernel of the homomorphism K ^ U ("H) is 
trivial. 

Consider the decomposition 

Sym{n) = 

The cone Ck{H) is related to the multiplicities through the following 
basic result, which is a particular case of Mumford theorem [ll](a proof of 
Mumford theorem, following closely Mumford argument, can be found in [8]). 

Proposition 24. We have m^(A) = 0 z/A ^ Ck{'H). Conversely, if \ is a 
dominant weight belonging to CKi'hC), there exists an integer k > 0 such that 
rn^{k\) is non zero. 

Thus the support of the function m^(A) is contained in the Kirwan poly¬ 
hedron CKifki) and its asymptotic support is exactly the cone Cxifl-L). 

Example 25. (The Cauchy formula, next) 
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Return to the example of the Cauchy formula. Let Ti = C”( 8 )C^ under the 
action of U (n) x U{N). The Kirwan cone has been determined in Example [TTl 
We see that the multiplicity function determined in Example [ 8 ] is supported 
exactly on the set A fl Cxi'H) (and with multiplicity 1). □ 

It is interesting to understand, for a given A G Ck{H) fl A>o, what is 
the smallest positive k such that rri^{kX) 7 ^ 0. We call this k the saturation 
factor. We will give one example, Example 1841 where k = 2. We would have 
liked to find larger fc, but we could not succeed. 

We now describe the nature of the function m^{\) on Ck{H). For this, 
we need to recall the definition of a periodic polynomial function. 

Consider a real vector space E with a lattice L. We can think of E = 
and L = Let C{L) be the space of functions on the lattice L. The 
restriction to L of a polynomial function on E will be called a polynomial 
function on L. Given an integer g, a function on L/qL will be called a periodic 
function on L of period q. If g = 1 we just have a constant function. 

Definition 26. A periodic polynomial function on L is a function on L 
which is a linear combination of products of polynomial functions with peri¬ 
odic functions. 

We also say that a periodic polynomial function p is a quasi polynomial 
function. 

The space of periodic polynomial functions is graded: if p(A) = j 
where the polynomials pj are homogeneous of degree k and the functions q 
periodic, we say that p is homogeneous of degree k. As for polynomials, we 
say that p is of degree /c if p is a sum of homogeneous terms of degree less or 
equal to k, and the term of degree k is non zero. If all functions Cj(A) are of 
period g, we say that p(A) is of period g. 

Example 27. 


is a periodic polynomial function of /c G Z and of degree 2. 

□ 

In other words, if p(A) is of period g, for any Aq, the function A —?■ 
p(Ao + gA) is a polynomial function on L. So we can represent a periodic 
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polynomial function as a family of polynomials indexed by L/qL. If q is very 
large, the above description is not efficient (the numbers of cosets being quite 
large). Nonetheless we will give an example of this description in Section 1631 

In practice, p will be naturally obtained as a sum of quasi-polynomial 
functions Pi,P 2 , ■ ■ ■ ,Pu of periods qi, q 2 ,..., qu- So p is of period q where q 
is the least common multiple of qi,q 2 , ■ ■ ■ ,qu- However, it is already more 
efficient to keep p as represented as YlVii fh® number of cosets needed to 
describe each pi being q^, and is usually much smaller that q. We thus 
will say that the set of periods of the quasi-polynomial function p is the set 
{?!, 92 , •••,?«}• 

An efficient way to represent periodic functions by step-polynomials is 
given in EH, i, as in the example we gave in the introduction. 

In the case of one variable, we can give the following characterization of 
quasi-polynomial functions p{k). If the function p{k) is quasi-polynomial, its 
generating series Taylor expansion for |t| < 1 of a rational 

function R{t) = where the a* are integers, and P{t) a polynomial 

in t of degree strictly less than The correspondence is as follows. 

Consider a quasi polynomial p{k) of period q, equal to 0 on all cosets except 
the coset f+ql^, with 0 < f < q. Write the polynomial function j —)■ p{f+qj) 
of degree R in terms of binomials: p{f -|- qj) = X]n=o Cn”) • Then 

oo R ^ 

(I — ■ 

j=0 n=0 ^ 

(Given a rational function R{t) = jp—algorithm (in polyno¬ 
mial time if the number s of factors is hxed) to compute explicitly the corre¬ 
sponding quasipolynomial function p{k) such that ~ -^(^) given 

in LattE integrate [5j.) 

As in the case of polynomial functions, to test if two quasi polynomials 
functions are equal, it is sufficient to test it on a sufficiently large subset. 

Lemma 28. Let pi,p 2 be two quasi polynomial functions. If there exists a 
open cone r, such that pi,p 2 agree on a translate s -|- r of r, then pi = P 2 - 

Recall that 

r = dimc('R) — |A^| — dimCx('H). 
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The following theorem can be considered as a quantum analogue of Duistermaat- 
Heckman theorem. (Theorem [2Uli 

Theorem 29. Consider the decomposition of the cone Ck{H) = CaCa, in 
the closed cones of polynomiality Ca, then, for each a, there exists a quasi 
polynomial function p^ „ of degree r on the lattice A such that z/ A G n A 


=PK,aW- 

This theorem is Meinrenken-Sjamaar theorem for the particular case of 
the action of K in the projective space associated to "H. In this work, the 
function m^{X) is related to the Kawasaki-Riemann-Roch number [27] (suit¬ 
ably dehned) of the reduced symplectic space 'Hred,\ = One 

may look in [15] for a different proof. 

Assume the cone CxiH) is solid, and let be the open subset of Ti 
consisting of elements with hnite stabilizers. Let Ca be a cone of polynomi¬ 
ality, choose a regular value G of the moment map, and let si, be the 
orders of the stabilizers for v G Then the set of periods of the 

quasi-polynomial function is contained in the set In particular, if 

the action of iL on "H has trivial generic stabilizer, the reduced spaces are 
smooth for regular values, and all the functions pf^- ^ are polynomials. 

Finding the set of periods is already quite non trivial. It is related to the 
computation of the Picard group of the reduced spaces EH- In the examples 
we study here, we compute a set containing the set of periods by brutal force, 
as we consider here relatively small values of the rank of K and the dimension 

oin. 

So, for any dominant weight A belonging to CxifH), the function k 
m^{kX) is of the form: m^(k\) = where Ci{k) are periodic 

functions of k. This formula is valid for any k > 0 (so co(0) = 1). The 
highest degree term for which this function is non zero will be called the 
degree of the quasi polynomial function m^{kX). We discuss this degree in 
the next subsection. 

If A is not in CK{T-i), the function k —)■ m^{k\) is just equal to 0, except 
for /c = 0, and conversely, if this function is not zero, then A G Ck{'H). 

An interesting particular case is when the center of 6 is one dimensional, y 
a weight of T vanishing on [t, t]nt, and such that y( J) = 1. So y indexes a one 
dimensional representation of K, and m^{kx) = The 

corresponding generating series is the Hilbert series HS(t) 

of the ring of invariant polynomials on "H under [K, K]. So, the degree of the 
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function k —?■ m^^kx) is the maximal number of algebraically independent 
invariants. There is a non trivial invariant P (different from a constant) if 
and only if 0 belongs to the Kirwan polytope of the projective space PilP)'. 
equivalently, if the line M>oX is an edge of the cone Ck{H)- This is one of 
the first instance of the [Q,R] =0 theorem, and this basic case follows from 
Mumford description [13] of the GIT quotient. The function k 'm^{kx) is 
a quasi polynomial on the full positive line fc > 0 (not only for k sufficiently 
large). It is in general difficult to decide if 0 belongs to the Kirwan polytope 
of the projective space T’('H), and even more so to determine the degree of 
the quasi-polynomial function k m^{kx)- 

The two following propositions are particular cases of Meinrenken-Sjamaar 
result. 

Proposition 30. Let (A,p, u) be a triple of dominant weights for U{d), be¬ 
longing to the Horn cone. The dilated Clebsh-Gordon coefficient k —)■ 
is a polynomial function of k for /c > 0. 

(If (A,/i, u) is not in the Horn cone, this function is identically 0, except 
if /c = 0). 

Indeed, we have already seen the quasi-polynomial nature of fc —)■ 
for /c > 0. Now it is easy to see that the generic stabilizer of the action 
of U{d) X U{d) X U{d) on H is the center of U{d) (embedded diagonally). 
Here is an explicit proof. Let Di,D 2 be two positive definite Hermitian 
matrices, with distinct eigenvalues. We assume Di diagonalizable in the basis 
ei, 02 ,..., Cd, while D 2 diagonalizable in a basis containing Consider 

A = D\^‘^, B = in 0 [((i) and {g,h,k) ^ U{d) x U{d) x U{d) such that 
gAk~^ = A, hBk~^ = B. We obtain AA* = Di = gDig~^, and BB* = D 2 = 
hD 2 h~^. Thus g commutes with Di, so is diagonalizable in the standard basis. 
So g commutes with A, and gAk~^ = A implies g = k. Similarly h = k, and 
h is diagonalizable in the basis diagonalizing D 2 . So we have g = h = k. 
But g being diagonalizable in the basis Cj, and in a basis containing '^ei, 
the equation g{ei 02 + ■■■ + Cd) = a{ei -f 62 + • • • + e^) implies that all 
eigenvalues of g are equal. 

A proof of Proposition [30] by more combinatorial methods is given in [T8] . 

We now consider the action of K = U{ni) x U{n 2 ) x ■■■ x U{ns) in 
0 • • • 0 C^^ 
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Proposition 31. Let A* G PKu{ni),>Q be polynomial dominant weights for 
U{ni) . We assume that (Ai,..., A^) belongs to the Kronecker cone CKifki)- 
The dilated Kronecker coefficient k g{kXi ,..., kXs) is a guasi-polynomial 
function of k for /c > 0. 

This result is asserted in |42j . 

(Remark that if (Ai,..., A^) is not in CxiTL), the dilated Kronecker co¬ 
efficient is identically 0, except if fc = 0). 

In the Kronecker case, we do not know the set of periods of the dilated 
Kronecker coefficient. 

Example 32. 

For the case ni = n 2 = = 3, we find quasi polynomials with set of 

periods {1, 2, 3,4} (see l6.2p . leading to polynomial behavior on cosets /-I-12Z. 

For the 4 qubits case ni = n 2 = = 2, we find quasi polynomials 

with set of periods {1,2, 3,4} fsee 16.2p . leading to polynomial behavior on 
cosets / -|- 6Z. 

For the 5 qubits case ni = n 2 = = 2, we hnd quasi 

polynomials with set of periods (1, 2, 3,4, 5} fsee l6.2p . leading to polynomial 
behavior on cosets / -|- 60Z. 

□ 

It would be desirable to describe as the number of integral points 

in a polytope. This result would imply directly the quasi polynomial behav¬ 
ior of m^{kfi). For Clebsh-Gordan coefficients, this is the hive polytope of 
Knutson-Tao [33]. The corresponding computation is implemented in im. 
For Kronecker coefficients, no general result is known. 

Let us end this section by recalling the behavior of multiplicities on the 
faces of 

Consider a face F of the cone Ck{H)- If F contains elements of CK{Ti) fl 
ft>Q, we then say that F is a regular face. If F is a regular face, let Tp be 
the torus with Lie algebra hn(F)^®''^. Let be the centralizer of Tp, with 
Lie algebra tp, and system of positive roots Let Ti{Tp) be the subspace 
of Ti stable by Tp. This is a Hamiltonian space for Kp. 

The following reduction formula holds. 

Theorem 33. Let F be a regular face of CKiTL)- Then for any X E F, we 
have 

m^(A) =m^f’^^(A). 
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Again, this theorem is a corollary of the [Q, i?] = 0 general theorem of 
Meinrenken-Sjamaar, (see also |1S]). 

Return to the example of the Horn cone. Let us consider a Horn triple 
(/, J, K) such that the equation (jl]) determines a facet of the Horn 

cone. Let F, the complement of /, J,K in [1,..., d]. Then the cen¬ 

tralizer of the element {Ej, Ej, Ek) is {U{I) x U{E)) x {U{J) x U{.E)) x 
{U{K) X U{K^)). Here we denoted by the unitary group acting on 

C/ = ©jg/Ccj. 

So when A,yU., z/ belong to the facet Ejj^k, the corresponding Clebsh- 
Gordan coefficient is the product of the Clebsch-Gordan coefficient 

relative to the group U (r) with the Glebsch-Gordan coefficient relative 

to the group U{d — r). So Theorem 15^ implies the factorization theorem of 
[28] . proved by combinatorial means. 

4.5 Comments on degrees 

We continue to assume that the moment cone Ck{'H) intersects the interior 
of the Weyl chamber and that the kernel of the homomorphism K — U{'H) 
is trivial. 

Recall the dehnition of 

r = dimc('H) - |A^| - dimGx('H). 

Gonsider Cq a cone for polynomiality and the two associated functions 
d^KaiP^a- The function is a polynomial function on it*, while is a 
quasi polynomial function on A. Then the two functions d^^ and p^a 
both of degree r. In the case where K is abelian, (partition functions), the 
term of p^ „ of highest degree r is polynomial. In general, this is not usually 
true. The term of p^ „ of highest degree r might be a quasi polynomial and 
not a polynomial. For example, r might be equal to 0, however the functions 
might be periodic function of period g > 1. 

The following example was communicated to us by M.Walter. 

Example 34. 

Let us consider U{2) and its torus T. Let it* with basis ei = [1,0], €2 = 
[0,1] and the lattice of weights A = {A = [Ai, A 2 ], Ai, A 2 G Z}. 

Gonsider K = f/(2)/{±l} and its torus T^- Then the lattice of weights 
is Ax = {A = [Ai, A 2 ]; Ai, A 2 G Z; Ai -|- A 2 even}. 
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Consider the 3-dimensional irreducible representation of K on H = 
(with trivial kernel), with weights [2, 0], [1,1], [0, 2], 

The Kirwan cone Ct{H) is M>oei © M>oe2. Let Ci be the closed cone 
generated by ei, ei + 62, and C2 generated by ei + €2, €2. 

The Duistermaat-Heckman function for Tk is given by 

+ ^2^2) = -^2 

on Cl, 

+ ^2^2) = -^1 

on C2. 

The multiplicity function for Tx on is (with Ai + A2 G 2Z) is 

+ -^262) = -A2 + ^ + 

on Cl, 

+ -^262) = -Ai + ^ + 

on C2. 

The (closed) Weyl chamber is Ci and we have CxiT^) = Ci- 
The Duistermaat-Heckman function is 

+ ^2^2) = 1 - 

The multiplicity function is 

+ '^2^2) = 

In particular, the highest degree term (degree 0) of the multiplicity function 
is not a polynomial. Only when Ai is even, we obtain that the highest degree 
term of (which is the constant function 1) coincide with the Duistermaat- 
Heckman polynomial. □ 

The degree r is equal to 0 if and only if multiplicities are bounded. This 
implies that all the fibers <h“^(.^) of the moment map are homogeneous spaces 
(the reduced fiber is a point). In particular, if the K multiplicities of the 
representation Sym{'H) are bounded, they take only the values 0 or 1. 
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We now consider a face F of Ck{'H). The restriction to F of the function 
is again a piecewise quasi-polynomial function. On F fl c^, the multi¬ 
plicity is given by the restriction of to F 0 Cq. The degree of Vk a\iin(F) 
(considered as a quasi-polynomial function A on the lattice lin{F) flA) might 
drop. If F is a regular face, we can compute the degree using Theorem [3^ 
Indeed 

Lemma 35. If F is a regular face, then the degree of the function 
restricted to F is the same on each cone Cq O F and is egual to rp = 
dime'H(Tp) - dim(F) - |A+|. 

If F is not a regular face, it is not easy to describe the degree of the 
restriction by such a simple formula. In particular if 0 is in the Kirwan 
polytope of P{'H), the degree on the corresponding edge of the Kirwan cone 
is difficult to compute. 

If F is a regular face, the degree of the multiplicity function restricted to 
any fl F depends only of the linear span of F. Let w be an element of the 
Weyl group. Although the Kirwan cone is not stable by the action of Wt, we 
might have two faces Fi and F 2 , so that their linear span is conjugated by w 
(but of course not the intersection of this linear span with Fi^('H)). We then 
obtain that the degree of the multiplicity function is the same on Fi, F 2 , see 
Example [821 

Consider now a line {/cA,/c = 0,l,2,...} contained in Ck{'H). It is clear 
that the degree of the quasi-polynomial function k —)■ m^{kX) is less or equal 
to r. The reduced space <h“^(FA)/F, which may be singular, is provided 
with a structure of projective variety. Thus we can compute in principle its 
dimension r^. More precisely, let us consider the open set of orbits of 
maximal dimension in F = $“^(FA). Then /K is an orbifold, and then 
rx = dim F°/ K. The following result is also a consequence of Meinrenken- 
Sjamaar. 

Theorem 36. The degree of the guasi-polynomial function k —)■ m^{kX) is 
egual to rx- 

However, it is usually difficult to compute explicitly rx using this geomet¬ 
ric theorem. We can do it when A is in the interior of the Weyl chamber. 

Lemma 37. Let X be in the interior of the Weyl chamber. Consider the min¬ 
imal face F of CKiTL) containing X. Then the degree of the guasi-polynomial 
function k —)■ m^{kX) is egual to rp = dimc'H(Tp) — dim(F) — |A^|. 
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Thus the behavior of the function is strikingly different from the be¬ 
havior of the Duistermaat-Heckman measure which typically vanishes identi¬ 
cally on the faces. So for a A which is not in the interior of the cone Ck{'H), 
the degree of the function dh^{tX) and of the function m^ikX) are usually 
different. 

A point A is called (weakly) stable, if the degree of A; —?■ m^{kX) is 0. In 
other words, the function m^{kX) is bounded on the ray with generator A. 
In this case the function m^{kX) takes only values 0 or 1. So this degree 
can be equal to 0 only on the boundary of "H, except in the case where 
multiplicities are bounded. Meinrenken-Sjamaar result (or in fact, here GIT 
theory) implies the following lemma. 

Lemma 38. A point X G Ck{H) is (weakly) stable if and only if ^~^{KX) 
is a K-orbit. 

A point A is called a stable point, if m^(A) = 1 , and if for any 7 G A>o, 
the function k -|- kX) is an increasing and bounded function of k. 

The following result, conjectured by Stembridge Ha, has been obtained 
recently by [IH]. 

Proposition 39. A weakly stable point is stable. 

L. Manivel [39] has determined some faces of the Kirwan cone Ck{'H) con¬ 
sisting of stable points in terms of compatible imbeddings, and has described 
the corresponding stabilized multiplicity hmfc_!.oo (7 + kX) in geometric 
terms. Explicit descriptions in the case of the Kronecker cone are given in 

Assume that the cone Ck{H) intersect t^g- If the inequations Xa > 0 
of the Kirwan cone Ck{'H) are known, the preceding discussion allows us 
to compute all stable points by elementary combinatorics, except on the 
boundary of the Weyl chamber, see Example [H2l 

4.6 The 3-qubits example 

Before explaining our method to obtain actual computations of the multiplic¬ 
ity function, let us give a complete simple example, HU, [In] . We followed 
the exposition of ra- 

We consider the action of 17(2) x f/(2) x f/(2) on H = ® ® C^. 
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The representation of U{2) x U (2) x U (2) in SymiTi) decomposes as 

Symin) = h, 0 

over polynomial irreducible representations A,/i, z/ of U{2) x U{2) x U{2). 
They are indexed as A = [Ai,A 2 ],/i = = [^ 1 ,^ 2 ] with Ai,/ij,r'j G 

Z, z = 1, 2, Ai > A 2 > 0, /zi > /Z 2 > 0, > 1^2 > 0. 

Though the dimension of a Cartan subalgebra of U{2) x 17(2) x 17(2) is 6, 
we can restrict our attention to the subset E C it* dehned by the equations 


E = {{a, 13, 7 ); «! + 02 = A + /32 = 7i + 72 }- 


Here a = [ai, 02 ], with ai, 02 real numbers, and ai > 02 , etc.. 

Indeed, considering the action of the center, it is easy that the Kirwan 
cone Ck{'H) is contained in E. More precisely, the two dimensional subgroup 



of 17(2) X 17(2) X U (2) is contained in its center and acts trivially on "H. So 
we might consider K = (17(2) x 17(2) x 17(2))/S' acting on "H, with Cartan 
subalgebra 1. Then it* is isomorphic to E. 

Higuchi-Sudbery-Szulc equations |23] (see Example [HD for the Kirwan 
cone in E are 


CKin) 


ai > Q!2 > 0, I 3 i > I 32 > 0, 71 > 72 > 0, 1 

0^2 < 72 + h , [32 <72 + ^2, 72 < 0^2 + /92 J ’ 


The image of pure states in "H is contained in the three dimensional affine 
space: 

El = {(a, (3, 7 ); «! + as = /3i + /32 = 7i + 72 = !}• 

The group S 3 of permutations of a,l3,'y acts on E. 

We parameterize E by by associating to [[ai,a 2 ], [/ 3 i,/ 92 ], [ 71 , 72 ]] the 
point (t, cq- 02 ,/di-/d 2 , 71 - 72 ) in with t = ai+ 02 = /Si+/S 2 = 71 + 72 - 
The Weyl chamber in = {(t,xi,Xs,X 3 )} is Xj > 0, and Bravy’s equations 
become 


{xi > 0, X2 > 0, X3 > 0, t + xi > X2 + X3, t + X2 > a;3 + xi, t + X3 > Xi + Xs}. 
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Figure 2: Kirwan polytope for three qubits 

As El is isomorphic to {(1,xi,0:2,xa)}, we can picture the Kirwan poly¬ 
tope A('H) in M'*' X R"*" x R"*". It is the polytope with 5 vertices ui, U2, ^3, "^4, ^5 
(see Figure [2]). 

For example, the point Vi := (0, 0, 0) corresponds to the triple Vmin in Ei, 
with 

= [[ 1 / 2 , 1 / 2 ], [ 1 / 2 , 1 / 2 ], [ 1 / 2 , 1 / 2 ]], 

and the point V 2 = (1,1,1) corresponds to the triple 

^rr^ax = [[l,0],[l,0],[l,0]]. 

. The other vertices n3,n4,n5 of A('H) are = [l,0,0],n4 = [0,l,0],n5 = 
[0,0,1], corresponding to the the triple [[1, 0], [1/2,1/2], [1/2,1/2]], and its 
permutations by S3. 

Thus the Kirwan polytope is an union of two tetrahedra glued over the 
triangle T with vertices [n3,n4,n5]. Consider the center c = (1/3,1/3,1/3) 
of this triangle, corresponding to the triple [[2/3,1/3], [2/3,1/3], [2/3,1/3]], 
and take the subdivision of the triangle T in 3 triangles, Ti,T2,T3. Then we 
have 6 cones of polynomiality: the 3 cones with basis the 3 tetrahedra which 
are the convex hulls of (Ti,Vmin) or the 3 cones with basis (Ti,Vmax)- 

Up to permutations by S3, we obtain two cones. Writing k = ai + a 2 = 

A + /52 = 7i + 72, 

Cl = {a,/9,7 I 0 < ai-Q;2 < mm(A-A,71-72) and ai-a2+A-A+7i-72 < iE\ 

and 

C2 = {a, /?, 7; 0 < ai-a2 < mm(A-A, 71-72); ai-a 2 +A-A+ 7 i -72 > 
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Figure [3] represents the cones C\ or intersected with the hyperplane Ei 
oiE. 



Figure 3: Cones Ci and C 2 inside the Kirwan polytope 

The multiplicities on the cone Cxi'H) have been described by Briand- 
Orellana-Rosas, m, [In] . On each cone g{X,^,h') is a periodic polynomial 
of degree 1 and period 2. Set fc = Ai + A2 = + 1^2 = hi + h2- Then 

u) = 

I i(Ai-A 2) + |(-l)"^+^^+"^ + i + 1 if (A,/i,z/) eCi 

1 Ai — I /ii + I A2 — ^ + A (—1)^^"'' ^ + I if (A, /i, z/) G C2 

In particular for example = [[!) 1]) [1) 1]) [1) 1]] is in Ci while Vmax = 

[[1,0], [1,0], [1,0]] is inC2. 

In this example, the multiplicity function takes values 0 or 1 on the bound¬ 
ary of CK{T-i). In particular, we see that g{[k, k], [k, /c], [/c, k]) = (—1)^|, 

so that g{\k, /c], [/c, fc], [fc, k]) is 1, or 0, if k is even or odd, and 

00 -| 

^ g{[k, k], [k, fc], [fc, k])t'^ = 
k=0 

as follows from the study of the Hilbert series of the ring of invariant poly¬ 
nomials on 0 (8) under the action of SL{2) x SL{2) x SL{2). 

We conclude this example with Figure 0] that shows the Duistermaat- 
Heckman measure for three qubits. The drawing is along the line [[1/2,1/2], [1/2,1/2], [1/2,1/2]] 
to [[1, 0], [1, 0], [1, 0]] between the bottom vertex Vmin and the top vertex Vmax 
of the Kirwan polytope np to the top vertex It grows linearly between Vmin 
and c, then decreases from c to Vmax- 
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Figure 4: Duistermaat-Heckman measure for three qubits: drawing along the 
hue [[1/2,1/2], [1/2,1/2], [1/2,1/2]] to [[1, 0], [1, 0], [1, 0]] 


Here is the discrete version of the multiplicity function. We compute 
g{\Qk + s, 6/c — s], [6/c + s, — s], [6fc + s, 6/c — s]) for s from 0 to 6. Here is 
the answer: 

[1,1, 3, 2, 2,1,1], forA: = l 

[1,1, 3, 3, 5,4,4, 3, 3, 2, 2,1,1], forfc = 2 

[1,1, 3, 3, 5, 5, 7, 6, 6, 5, 5,4,4, 3, 3, 2, 2,1,1], for A: = 3 

It shows clearly the quasi polynomial behavior of the multiplicity on cones 
of polynomiality. 

Our program compute symbolically the multiplicity function on the do¬ 
mains of polynomiality. 

4.7 Multiplicities and Partition functions 

We are now going to explicit the quasi polynomial functions determining the 
Duistermat-Heckman measure for the space "H, and the multiplicity function 
for Sym{'H) whose existence is stated in Theorem [20] and Theorem [29l 
For computing multiplicities, the connection is made trough partition 
functions, whose computation is achieved (using techniques developed in 
0 ). via iterated residue of rational functions with poles on arrangement 
of hyperplanes. We give a simple example, explaining the philosophy of the 
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method, and we compare it with the algorithm of Guoce Xin [SI] in Snb- 
section 14.9.91 The reader may want first to read this subsection. We state 
formulae Theorem [51] and Theorem 1561 for the Duistermat-Heckman measure 
and Theorem l53l and Theorem l5^ for the multiplicity. Even if we could follow 
the same pattern, we will take a different approach for computing Kronecker 
coefficients. We will compute them as a byproduct of the computation of 
branching coefficients, as explained in detail in Section 0 

To start with, when K = T is abelian, computation of multiplicities 
is the same problem than computing the number of integral points in a 
(rational) polytope, while computation of Duistermaat-Heckman function 
is the same problem than computing the volume of a polytope. In turn 
volumes and number of integral points in polytopes can be computed, as we 
just said, as iterated residue of rational functions with poles on arrangement 
of hyperplanes. The precise statements are collected in Theorem [5T] and 
Theorem (53] 

We write 

Sym{'H) = ©^gy,m^(/i)e^. 

Let T C it* be the list of the weights for the action of T on "H counted 
with multiplicities: 

= ['01,'02,---,'0iv]- 

We assume that the cone Cone{'$) generated by T is a pointed cone: 
Gone(T) n —Gone(T) = {0}, and that the lattice of weights A is generated 
by T. 

Let Pij, be the function on A that computes the number of ways we 
can write /i G A as ^ Xi'ijji with Xi nonnegative integers. The function 
is called the Kostant partition function (with respect to T). It is 
thus immediate to see that 


The moment map '■'H ^ it* is given by 

N 

^xi^ZaCa) = 

a=l 

Here we have used an orthonormal basis ei,e2, ... ,eN of Ti where T acts 
diagonally with weights [-01,...,-^tv]. So the cone Ct{'H) is just the cone 
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Cone{^) generated by the list \1/ of weights. Assume for simplicity that \1/ 
generates it*. It is thus a cone with non empty interior in it*. 

Let y E it*. Dehne the polytope 


N 

nviy(|/) = {[xi,...,Xn] e R^,Xi > 0, Xgipa = y}. 

a=l 

The Duistermaat-Heckman function dhj;{y) is the volume of the polytope 
n^(?/) and (/i G A) is the number of integral points in the polytope 

nv|/(p,). 

Dehne to be the open subset of the cone Criy-), where we re¬ 

moved from CTiT-L) the union of the boundaries of the cones generated by 
all subsets of T. Let us write = UYa where Ya are the connected 

components of Dehne Ca = Ya, so we obtain a cone decomposition 

CxiH) = UaCa- 

Here all cones Ca are closed and have non empty interiors. 

The following result explicit Theorem [20] in this setting and the well 
known behavior of the partition function . 

Proposition 40. • The function dhfp{y) is given by a homogeneous poly¬ 

nomial function d^ on Ca. The degree of d^ is egual to |T| — dimt. 

• There exists a guasi polynomial function Pa on A such that 

=pt{T) 

on the closed cone Cq. The degree of Pa is egual to |\k| — dimt and its 
highest degree term is egual to df^. 

• In particular m^(p) = Paip) on the closed cone c^. 

Let A be a compact Lie group acting on Ti by unitary transformations, 
and let us assume that the action of its maximal Cartan subgroup T has 
hnite multiplicity. In other words, the weights of T in "H span a pointed 
cone. Then 

Sym{H) = 
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and we have, for A G A>o, 

X] e(M;)Pvi,(A + pt-M;p*). 

uiEWj 

The cone CK{T-i) is contained in CT{T-i) H it>o hut usually smaller. Let 
us give some simple examples. 

Example 41. 

Let us consider the standard action of U{2) on with weights ei = 
[1, 0], €2 = [0,1]. The moment cone for T is thus M>oei © M>oe2. The Weyl 
chamber is t>o = {^iCi +^262; ^1 > ^2}- We have = zz*, if we consider 

2; G as a map —?■ C. Thus ^k{z) is a Hermitian non negative matrix 
of rank 1. We thus have Ax('H) = M>oei. In particular, the cone Ct('H) is 
solid in t*, while CxiH) is not solid. □ 

Even when both Ct{T-L) and Ck{T-L) are solid, they can be different. 
Example 42. Ct{'H) for the 3-qubits. 

Return for example to the case of 3 qubits 14.61 and we keep the same 
notations. The weights of the action of T in "H = C2 © C2 © C2 on the 
weight vectors e± © e± © e± with = Ce+ © Ce_ are [[1, 0], |1, 0], [1, 0]] and 
its permutation by the Weyl group S2 x S2 x S2 of K. In the parameters 
the weights are [1,61,62,63] with 6j G {1,-1}. So we see that C'r('H) 
is the cone over the cube C with vertices [61,62,63] with 6* G { — 1,1}. The 
intersection of CTiT-L) with the positive Weyl chamber is the cone over C 
intersected with the positive quadrant. This is again a cube with vertices 
with izi G {0,1}. 

Figure 0 shows the Kirwan cone CxifH) inside CTifH) nit>o. We see that 
the points (1,1, 0) (and permutations) are vertices of CTifH) nR>0) but they 
do not belong to CKifhi). The Kirwan polytope Ck{H) is obtained by re¬ 
moving from CTi^H) a ’’neighborhood” of these vertices, namely the simplices 
with vertices (1,1, 0), (1,1,1), (1, 0, 0), (0,1, 0) (and its permutations). 

□ 

Assume (to simplify) that K acts on T-L with generic hnite stabilizer. 
Then CKi'H) is a cone in R>q with non empty interior and it follows from 
the formula for Duistermaat-Heckman measure that the Kirwan cone Ck{'H) 
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Figure 5: iF-Kirwan polytope for three qubits inside T-Kirwan polytope 

is the union of the solid cones Canil>Q contained in it. As seen in the preceding 
example, it is not easy to understand which ones have this property. 

We denote by d_af{y) = ^f{y — ta)\t=o the derivative of a function 
/ on ii* in the direction —a. Let (n„>o d-o) the product of derivatives 
with respect to all negative roots in Af. Consider the cone decomposition 
Cr('H) = UaCa. 

Theorem 43. Assume CKi'H) is a solid cone. 

Let Ct{T-L) = UaCa. Consider Ca such that Ca n it>o C Ck{LL). 

• On Ca n it>Q we have 


o>0 


• Let Pa he the quasi polynomial function on K such thatV^^,{|a) =pt{p) 
on Ca- Then, on Ca n A>o, we have 

^{w)pl{\ + pi-wpt). 

uiSWj 


Remark 44. 

We see that the degree of the polynomial function dh^, and on each 
Ca is at most equal to the degree of the corresponding functions for T minus 
the number |A^| of positive roots. Our algorithm will indeed be easier for 
K than for T. 

□ 
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Proof. Consider first the open snbset tq of fl it>Q, at distance Q of its 
bonndary. We assnme Q greater than the norm of \\pi — wpi\\ for all w G Wf. 
If we use the formula 

= X] e{w)Vm{\ + Pi- wpk) 

and A is in tq, all points A + pt — wpk are in Ca and we may replace by 
the function We thus obtain on tq the expression 

= '^^{w)p^{X + pi-wpk). 

wGWi 

However the function A —)■ + Pi — wpi) on the right side 

of the above equation is a quasi polynomial function of A on Ca n (H>o)- 
By Meinrenken-Sjamaar, Theorem [29l we know that m^(A) coincide with a 
quasi polynomial function on the closed cone Cafl A>o, under our assumption 
that Ca n it>Q C CxiPL). 

Using Lemma |28l we obtain the equality on Ca n A>o. □ 

4.8 Back to degrees on faces 

Assume that the cone CxiPL) is solid. We rewrite the formulae for the de¬ 
grees rp on faces of the multiplicities, using T. Let F be a facet of CxfH) 
determined by an equation X = 0 intersecting (that is X is not propor¬ 
tional to a fundamental coroot Ha). Then we compute 'H(X), the subspace 
of PL annihilated by X. Its complex dimension is the cardinal of Tq, where 
Tq is the sublist of T composed of the ipj such that fjjiX) = 0. The system 
Aq of positive roots for t{X) is the sublist of Aj*" composed of the roots aj 
such that djiX) = 0. So we obtain the formula 

Lemma 45. On the facet F, the degree of the multiplicity function is |To| — 

I Aq I — (dimt — 1). 

4.9 Tools for computations 

From now on, we focus on the computation of 
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4.9.1 Residue calculations: the philosophy 

Consider the case of n dimensional torus T = {(mi, M2,..., Wn); l^il = 1} 
acting onT-L = . We denote by m —)■ = Hifhe character of T indexed 

by A G A = Z"^. Assume that all weights ipi of the action of T on have 
non negative coordinates. Almost by dehnition, the computation of 
is the computation of the Fourier series of the function S{u) = 

I Mil < 1. So we obtain, choosing e, small positive real numbers. 


N 


771: 


H 


(A) = 


M 


-A 


'\uj\=ej 


n 1— u'i’j n 




duj 
2i7ruj 


(5) 


We choose an order of integration, and the residue theorem in one variable 
allows to integrate the variable Mi, but we obtain a priori A^ — 1 new integrals 
in n — 1 variables, by considering all poles in mi of the functions (1 — m'^-’). 
Not all of these poles are in the interior of the circle |mi| = ei, then we have 
to keep track of the branchings. At the end we obtain a very large number 
of products of one dimensional residues computations. 

The essence of Jeffrey-Kirwan residue is to compute a priori the paths 
which contributes. In fact, it is easy to see that the function 
be written as a sum of functions fa, where each fa depend of n independent 
variables. Here a varies in a large set T. For example 


1 

(1 - Mi)(l - MiM2)(1 - M 2 ) 


1 

(1 -MiM 2)2(1 - Ml) 


+ 


_ U2 _ 

(1 - MiM 2)^(1 - M2 ) 


and we can separately use the independent variables (mi,miM 2) or (m 2 ,miM 2). 
In the residue computation at A, only a certain T(A) of these functions /„ 
will contribute. So we compute a priori what is the subset T(A) and we are 
reduced to product of residues in one variables. This is the object of the 
computation of adapted Orlik Solomon bases. We give now the technical 
details of the resulting formula. 


4.9.2 Definitions of regular elements and topes 

Consider a hnite set F of non zero vectors in a real vector space E. 

Definition 46. We say that a hyperplane H <Z E is F-admissible if H 
is generated by elements of F. We denote by A{F) the set of admissible 
hyperplanes. 
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Assume that F generates a lattice L in E. Let a be a subset of F, 
generating F as a vector space. We say that a is a basis of F. We denote by 
da the smallest integer such that d^L is contained in the lattice generated by 
cr. 

Definition 47. Consider the set periods{F) = {d^} where a runs over 
the basis of F. Let q{F) is the least common multiple of the integers in 
periods{F) We say that q{F) is the index of the finite set F (with respect to 
L). 

Definition 48. Consider a finite set F of hyperplanes in a real vector space 
E. We say that E E is F regular, if f is not on any hyperplane belonging 
to F. 

Let EregT be the set of F-regular elements. A connected component r of 
Eregj^ wHl bc Called a tope (with respect to F). If f is F-regular, we denote 
by r(^) the unique tope containing 

Return to the situation where 4^ C it* is the set of weights of T in IL. 
When F is the family .4,(4/) of 4/-admissible hyperplanes, we also say that f 
is 4/-regular instead of F-regular. A tope for the family 4.(4/) will be called 
a 4/-tope. 

4.9.3 Iterated residues and Orlik-Solomon bases 

A list = ['011) '012, • • •) i^ir] 4 elements of 4/ will be called a ordered basis 
if the elements '01^, form a basis of it* . Let ^(4^) be the set of ordered bases. 
An ordered base = [0ij, ..., 0i^] is an Orlik Solomon base, OS in short, 

if for each 1 < / < r, there is no j < ii such that the elements , ipi^ 

are linearly dependent. Denote by 05(4/) the set of OS bases. 

Definition 49. If T is a ^-tope, we denote by 05(4^, r) = {a G 05(4^), r C 
Oone(cr). Here Cone{a) = 

The set 05(4/, r) is called the set of OS adapted bases to r. 

We will show in Section 14.9.61 how to compute adapted bases. For an 
account of the theory cf. [2]. 

We now dehne the notion of iterated residue. If / is a meromorphic 
function in one variable z, consider its Laurent series at z = 0. The 

coefficient of z~^ is denoted by Res2=o /• 
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Let 0^ = [«!, 02 ,..., Or] be an ordered basis of it*. For z G tc, let 
Zj = (z, Oj), and express a meromorphic function f{z) on tc with poles on a 
union of hyperplanes as a function f{z) = f{zi, Z 2 , ■ ■ ■, Zr), (in particular / 
may have poles on zj = 0). 

Consider the associated iterated residue functional dehned by: 
Res-^{f{z)) := Res^i=o(Res^2=o • ■ ■ (Res^^=o f{zi, Z 2 ,..., Zr)) ■■■). (6) 


4.9.4 Multiplicity and Duistermaat-Heckman function for the 
torus T 

Fix a cone of polynomiality Cq, then is the union of the closures of the 
topes r contained in Cq. 

Topes were also called small chambers in [1]. 

We might need several topes (small chambers) so that the union of their 
closures is Ca- Figure [H] shows an example for the chambers complex versus 
the topes complex for the cone associated to a root system of type ^3. 


e-e, 



e-e. 



Figure 6: 8 topes (left) versus 7 chambers (right) 


Let G R*. Dehne the following function of 2; G tc- 

1 






The function z —)■ z) has poles on the union of the hyperplanes 

{'ijjyz) = 0. An iterated residue of the function z —)■ s^(^,z) depends of ^ 
through the Taylor series of the function at 2; = 0. So if we consider 
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^ as a variable, we obtain as a result of performing an iterated residue a 
polynomial function of 

Definition 50. Let t he a ^-tope. Define 

dri0= ReSaST{^,z). 

o-G05('I',t) 

Thus df is a polynomial function of 

Theorem 51. Let r C it* be a ^-tope, and ^ G r. If t is contained in 
Cone{^), 


Let p G A, and for z ^ ic, define : 


St ( p -, z) 


— pifj'D _ 


Let r be the dual lattice to A. Let q := g(T) be the index of T. So if a 
is a basis in it* consisting of elements of T then qA C 

If 7 G r, and we apply an iterated residue to the function z ^ z + 

^7), we obtain a quasi-polynomial function of p. Indeed, it depends of p 

through the Taylor series at z = 0 of and 

is a periodic function of pi of period q. 


Definition 52. Let t be a ^>-tope. Define 

PriL)= Y1 Y1 Res„S^{p.,z+—-f). 

cre05(^,r) 7er/gr ^ 


Theorem 53. (Szenes-Vergne) Let r C it* he a "h-tope and p, ^ r. If r is 
contained in Coneifi!), then for any p G r fl A, 

=pf(A)- 


Remark 54. 
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We denote by m —)■ the character of Tc associated to G A. Consider 

the function F{u) = clear that is the Fourier coeffi¬ 

cient of the function F{u), expanded in the domain \u'^\ < 1, for all V’ ^ 
Thus rnip^jj) = Here r is a small positive number. 

Using the change of variables Ui = (z* being the coordinates on it 

associated to the basis a), an iterated residue Res^s^{^, z) can be computed 
as an iterated residue at m = 1 G Tr of the function tt——=vt. 

Similarly the other terms associated to 7 G F/gF can be computed as 
iterated residues at some elements of hnite order q in Tc- Thus Szenes- 
Vergne theorem is a multi-dimensional residue theorem. A detailed example 
will be given in 14.9.91 □ 


4.9.5 Multiplicity and Duistermaat-Heckman function for K 

We now want to compute the functions dh^{^) and 
Let ^ G ft*. Dehne the following function of z G tc: 




npevi/(V',^) 


Definition 55. Let t be a '^-tope. Define 


<7e05('i',T) 

Theorem 56. Let r be a '^-tope and let ^ G rnft>o. //rflt^o is contained 
in the Kirwan cone CxiLL), then for any f Er H t>o, 

= Dpo 


Remark 57. 

The function Dfi is easier to compute that df as the orders of poles in 
of the function s|-(^, z) are smaller that those of sf (^, z). 

□ 

In particular, if ^ G ft^g; T-regular, we can check if ^ is in CxiR), 

by computing Dfi where r is the tope containing 
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We can now state the result to compute m]^(A). 
Let A G A>o, and for z G tc define: 


S^iX.z) 




n„,A+(i 


Definition 58. Let t be a tope. Define 

0,0 /jr 

^r(A) = Res^SK{\,z + -7). 

creOS{^,T) ■y&/qr ^ 


Theorem 59. Let t be a tope and let A G A>o. Then 

1. If X & T, then m^{X) = Pr{X). 

2. if X Et and r fi it>o C CxiTL) then mf^{X) = Pt{X). 

Proof. Consider the quasi-polynomial function, defj52) 

pf(T)= Y1 Y1 Res^S^{fr,z+—-f). 

crG05('I',r) 7 Gr/gr ^ 

The function Vi^, is given on r by the quasi polynomial formula, Theorem |53] 

=pf(/i). 

We now use Theorem 03] which asserts that is given on the closure of r 
by (/i + pt- wpk). This is 

X] e(M;) Y Res^'^T(h-'«^(Pj)+P{,^ +—7) = 

w&Wt aeOS{^,T) 7Gr/gr ^ 


^ ^ Res^^e(M;)e<-“'(^b+P«,-+i^7> 

(TGC><s('f ,t) 7gr/i}r w 


g(M,z+2^7> 



Y Y —7). 

f7G05(^,r) 7Gr/gr ^ 


The last equality follows from the denominator formula e{w)e^‘ = 

no,>o(l ~ e“)- Thus m^{p) = Pfi{p) when p is in the closure of r. □ 
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4.9.6 OS bases 


Let f be a \[^-regular vector. Let n = dimt. We recall here the algorith¬ 
mic method based on the method of Maximal Nested Sets (MNS) of De 
Concini-Procesi, [16], and developed by Baldoni-Beck-Cochet-Vergne in |2] 
to compute 05(4^, r(n)). Recall that C>iS(\l', r(n)) depends only of the tope 
t{v) where v belongs. We hrst need to compute the set A of admissible 
hyperplanes generated by elements of \1/ (we do not have a good method to 
perform this step). For each hyperplane, we choose an element FT G t, so 
that our hyperplane is H-^. 

Anyway, assume that we have determined this set A of admissible hyper¬ 
planes. Then we order T. We compute recursively a set of OS basis for T 
by the following method. 

• We hrst choose "00 the lowest element of T, and an hyperplane Hi 
where 'ipo does not belong: {Hi,'ipo) ^ 0. Let Ti = T fl H^ and 
Ai = e A, {H, 'iJjq) = 0 }, that is Ai is the set of hyperplanes in A 
that contain 'i/jq. 

• We choose '0i the lowest element in Ti and H 2 in Ai , such that 'ipi does 

not belong to H2, and such that 0 is spanned by T 0 0 Ff^. 

In practice, we verify this condition only at the end, but we verify that 
the set T 0 fl iF^ has at least n — 2 elements. 

We can now dehne T2 = fl iF^ O T and A 2 = {FT-*- G Ai, (FT, -^1) = 
0}, that is A2 is the set of hyperplanes in A that contain '0o, V'l- 

• The inductive step is the following: 

suppose we have constructed a list [ao,ai,... ,as] C T and a list 
[Ffi, H2, , Hg^i] of hyperplanes with the property that 

— Qfj is the lowest element of T fl n ■ ■ ■ fl H^. 

- aj ^ Hj^i and G Aj where Aj = {H^ G A, {H, 'ifjt) = 0,t = 

Thus [ao, ai,..., a^-i] G and aj ^ Hj+i. 

- T n Hj^ n H^ n ■ ■ ■ n Hj- n Hj:^^ has > n — s — 1 elements. 

• The next step is as follows: 

- we choose the lowest element in TflFfj^niF^n- • ■nFf^nFf^;^ 
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- We choose an hyperplane i?s+2 containing [oq, «!,..., a^], i.e Hs +2 ^ 
hut not Os+i and such that \1/ fl fl fl • ■ ■ fl Hj- fl 
has > n — s — 2 elements. 

• we then continue with [Hi, H 2 , ^^5+2] and [oq, ai,..., Os+i, as+2] 

• At the end, we obtain n elements [oq, ai,..., an-i]. 

It may occur that we cannot go to the end. This means that at this 
step, the space was not generated by its intersection 

with T. However, if we arrive to the end, we obtain a basis, and by 
pT] we obtain a set of OS basis for T. 

We can rehne this method to obtain directly the set of OS adapted basis 
for a regular vector v. The algorithm we just explained is modihed easily. At 
each step when we look for the hyperplane Hg containing [oq, «i,..., a<i-2], 
but not Os-i, we also impose on Hg the condition that the vectors v and 
ag-i lie on the same side of the hyperplane Hg. In this way, we obtain 
OS{'^,t{v)). If V is not in the cone 0(11/), then the algorithm returns the 
empty set for C>iS(T, r(n)). 

Remark that if we know the equations of the cone 0(\k), then we could 
hrst check if v is in the cone, before computing (9iS(T,r(n)) to shorten the 
procedure. However, we do not use this preliminary step. 

4.9.7 Dilated coefficients, Hilbert series 

Once we know (95(T,r) the calculation is not more difficult to do with 
symbolic variable p, and we obtain the periodic polynomial which coincides 
with m^(p) on r n it>o (and the closure if r fl H>o is contained in Ck{H)) . 

In particular, to compute the function m'^{kij) on the line Mp as a pe¬ 
riodic polynomial function of k is not more difficult than to compute the 
numerical value 

A particular interesting example is the computation of Hilbert series. 
Assume that t = ^©[t, t], where the center 3 = MJ of i acts by the homothety. 
Consider y G A such that y(i J) = 1 and y = 0 on z(tn[C t]). Then m^(/cy) = 
So the series R{t) = is the Hilbert series 

of the ring of invariants polynomials under the action of [Kc,Kc]. It is of 
the form 
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The degree of the quasi-polynomial function k {'H){kx) as well as 

its set of periods gives some information on the number of factors and the 
coefficients aj in this Hilbert series. 

4.9.8 Advantages and Difficulties of the method 

Since one of the objective of this paper is to describe an efficient algorithm 
to compute multiplicities, at least for low dimension, let’s look at the weak 
and good points in implementing the above formulae. 

• The computation of the set of admissible hyperplanes for a system T. 

For example, if we consider the space of 6 qubits, a brute force computation 
(taking any subset with 6 elements of 2® elements would lead to 7624512 
computations). We do not know an efficient algorithm for computing the 
set A(T). Moreover, our system T has symmetries coming from the Weyl 
group action. So it would be good to have an algorithm computing A(T) 
up to Weyl group action. When T is a root system of rank r, this set, up 
to Weyl group action, is just the set = l,...,r} where hi are the 

fundamental coweights. In the Kronecker examples, we do not know the set 
A(4/), although the system T is quite simple, see Formula [H] inside Example 

Eni 

• The residues calculation for the function z). We do this by power 

series expansion, and it leads to multiply polynomials of larger and larger 
degree. 

• The computation of the set OiS(T,r). Once this set is computed, 
residues contributions are independent of each other. So an advantage of the 
iterated residue method is that each individual residue computation is easy 
and does not use much memory. 

• The index q of the set T. The fact is that for arbitrary systems T 
with iV-vectors, this index can be large. So Szenes-Vergne formula does not 
provide a polynomial time algorithm (the dimension of H being fixed). 

For example, in the case of the knapsack, (Example [7D, when the Aj are 
relatively prime, the index g is Ai • • • A„. Thus, if q is large, the function 

>S'|(/i, ^ + —7) 

should be computed in polynomial time using Barvinok determination of 
generating functions of cones [6]. This is the method we used in [5] to give 
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an efficient computation when n is fixed. However, for our computation 
of Kronecker examples, we just used the summation over L*/qL*, as q was 
relatively small. 

Let us explain at this point the method followed by Christand-Doran- 
Walter [13]. Their method (of polynomial complexity, when the dimension 
of H is fixed) computes the numeric multiplicity by the following method. As 
we said, the multiplicity function rnip{\) for T is the number of points in the 
polytope n,j,(A). Thus [13] uses the ’’Barvinok algorithm”, as implemented 
in [51], to compute m^(A). This computation would also be possible to do 
using the Latte software package [5]. So the above difficulty of the possibly 
very large index q is resolved (in polynomial time) by the power of Barvinok 
signed decompositions. Then they use the relation between and to 
compute the multiplicity under K. Remark that this last computation can 
be done on parallel computing. However the number of elements w becomes 
very large. 

We close this subsection with some questions. 

Questions 

Here are the problems that we would like to have some partial answers in 
interesting examples: 

• Can we understand the image of the moment map, and describe the 
Kirwan polyhedron by inequalities? 

• Can we compute the Duistermaat-Heckman measure? 

• Can we compute the dilated multiplicity k —)■ m^{kX) = 

or at least say something on the periodicity of the coefficients and the 
degree of the function k —)■ m^{kX) ? In particular when A is a one 
dimensional representation of K, or other interesting A. 

4.9.9 A very simple example 

Return to Example 14.51 We close this part in computing in two ways the 
multiplicity m^^(A) in this very simple example. 

• The straightforward method of computing the expansion of a Fourier 
series. The one dimensional residue formula is used repeatedly to replace the 
large poles at u = 0 in poles on |m| = 1. The branchings appearing in this 
straightforward becomes soon quite complicated, if the rank of K as well as 
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the dimension of Ti increases. This is the method followed by [M], with an 
efficient way to keep track of the branchings. 

• The Jeffrey-Kirwan residue method. It gives a (possibly large) number 
of residues at 2; = 0, or equivalently, in exponential coordinates u = exp z, at 
M = 1, or some hnite order elements in T. The computation for each iterated 
residue are independent of each other. 

• The straightforward method. 

We consider the torus T of t/(2) parameterized as {(ui, M2 ); |mi| = 1; |w 2| = 
1} and the action of T on "H with weights [2,0], [1,1], [0,2]. The character 
of the representation of T in Sym{'H) is obtained by computing the Fourier 
expansion of 

F{Ui,U2) 

for mi,M 2 of modulus strictly less then than 1. Let A = Aiei + A2e2, with 
Ai, A2 two non negative integers. Thus, taking for example ri = r2 = |, we 
obtain the straightforward formula 


m^(A) = 



u 




dui du2 


I„l|=r-i>2|=r2 (1 -UiU2){l - «!) 2i7rMi 2i7rM2 ' 


Let us fix Ml, and integrate on |m 2| = 'r2. As A2 > 0, the poles of the 
integrand on |m 2| <r 2 is the pole at M2 = 0. The other poles are M2 = 1, Mi = 
—1,M2 = 1/mi, and there are no poles at 00. Using the one dimensional 
residue theorem, we obtain (A) = A + B + C with 

^ If dui 

2 (1 - wi)(l - ul) 2i7rui ’ 


B = -tDl f dm ^ 

2 (1 - Mi)(l - M?) 2i7rMi’ 

r dMi 

i|m|=ri (I - ujy 2iTTUi' 

We now integrate in mi. However, at this step, the computation is different 
if A2 > Ai or not. Indeed in the case where A2 > Ai, the integral C will be 
equal to 0, as the integrand has no poles inside |mi| < So let us assume 
this is the case. We then pursue in computing similarly the integrals A,B by 
the one dimensional residue formula, and we obtain by summing the residues 
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at Ui = l,Mi = —1 of the integrands expressions in A,5, a sum of 4 terms 
adding up to 

^(1 + + (1 + + ^((- 1 )"^ + (- 1 )"^), 

an expression vanishing if Ai + A2 is odd, as it should, and coinciding with 
the given expression 


on the cone of polynomiality C2 

• Let us now employ Theorem ESI using Jeffrey Kirwan residues. 

We consider the same case A2 > Ai. 

Consider the tope T 2 = {A; Ai > A2 > 0}, the interior of C2. If we compute 
OS'(\l/, T 2 ) for the order T = [[2, 0], [1,1], [0, 2]], there is only one adapted basis 
a = [[2,0], [0,2]]. The index q is equal to 2 (in the lattice generated by 
T). The function St{\,z) is equal to 

gXlZi+\2Z2 

(1 - e-2^i)(l - e-2^2)(i _ e-^ 1 -^ 2 ) ■ 

A representative of Tx/STx is G = [—1/2,0], and we obtain 

St{X, z + i'K'y) = Si + S2 

7 GrK/ 2 rK 


with 

„\lZi+\2Z2 

c ^ _r_ 

^ ^XlZl+\2Z2 

*^2 = (- 1 ) ' (1 _ e- 2 -i)(l - e-2-2)(i + e--i--2)' 

The iterated residue computation ReSz^=QReSz 2 =oSi is straightforward 
and we obtain |Ai + | from Si, and |(—1)^^ from the second term 5*2. 

It is also worthwhile to remark in this example that the iterated residue 
computation depends of the order. The reverse order Resz 2 =oR€-szi=oSi would 
have given the formula for the tope ti = {(Ai, A2), Ai > A2 > 0}. 
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5 Branching Rules 

In this section, we will write a qnasi polynomial formnla for the branching 
coefficients. This will allow, in combination with the Canchy formula, to 
express the Kronecker coefficients in a different way that has at least the 
advantage of reducing the number of variables in the setting. Section |6l 

5.1 The Branching Cone 

Consider a pair K <Z G oi two compact connected Lie groups, with Lie 
algebras t, g respectively. Let tt : g* —)■ P be the projection. 

We consider the following action of G x iL on G : {g,k) ■ a = gak~^. The 
manifold T*G (the cotangent bundle of G) is a G x iC Hamiltonian manifold. 
The geometric quantization of T*G ”is” the space L‘^{G). This statement 
can be justihed, but we will not do it here. 

Let us dehne V = R{G) to be the subspace of G°°(G) generated by the 
coefficients {gui, U 2 ) of hnite dimensional representations of G (by Peter-Weyl 
theorem, the space L‘^{G) is the Hilbert completion of V). 

Assume G,K connected, and let Tg, Tk be maximal tori of G,iC. We 
may assume, and we do so, that Tk C Tq- We choose Cartan subalgebras 
tg, h, Weyl chambers >q, and we denote the corresponding cones of 

dominant weights by Ag,>o, Rk,>o- We denote by ^ >q the sum H* >q©H^>q 
of the closed positive Weyl chambers relatives to G, JL, and by H* (>0 i^s 
interior. We may also choose compatible root systems on K, G: If A is 
dominant for G, then the restriction of A to is dominant. 

For A G Ag >0 (resp. n G Ax,>o), denote by (resp. V^) the irreducible 
representation of G (resp. K) of highest weight A (resp. p). 

Under the action of G x G, 

So, under the action of G x iP, 

V = ©A,/.mG,i^(A, p)Uf (g) 

(A varies in Ag,>o, and p in Ak,>o) where g) is the multiplicity of the 

representation /i in the restriction of Up to iP; it is also called the branching 
coefficient. 
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Let us write coordinates for T*G. We identify the tangent bundle TG 
to G X Q through the left translations: to {g,X) G G x g we associate 
g\t=o e TgG. Thus, T*G is identihed to G x g*, and the moment map 
relative to the G x iL-action is the map ■ T*G —)■ g* © 6* dehned by 

{g, —)■ —7r{g~^ '0)- order to relate the moment map to the branching 

coefficient we use the slightly modihed map: <h : T*G — )■ *g*©zt* 

given by 

^ -O). 

Here g & G, and ^ G g*. 

Let 

GG,K{T*G) = ^T*G)mtl,^^, 
be the Kirwan cone associated to <h: 

CaAT’G) = e !t;>„ x ,, e 7r(G ■ 0} • 

The set Gg,k{T*G) is a polyhedral cone in f >q, and is related to the 
branching coefficients through the following basic result. 

Proposition 60. VTe have mG,i<:(A,/i) = 0 if {X, fi) ^ Gg,k{T*G) . 

Conversely, if {X, fi) is a pair of dominant weights contained in Gg,k(T*G), 
there exists an integer k > 0 such that mG,K{kX, kfi) is non zero. 

Thus the support of the function mG,K{X, fi) is contained in the Kir¬ 
wan polyhedron Gg k(T*G) and its asymptotic support is exactly the cone 
CG,KiT*G). 

Remark that if G = iL, the cone Gg,k(T*G) is just the diagonal 0 ) ^ ^ 

it* >q} in it* g >Q. However, we assume from now on that no nonzero ideal 
of t is an ideal of g (this condition excludes the preceding case). It implies 
the following result (Duflo, private communication). 

Lemma 61. The polytope Gg,k(T*G) is solid. 

Let TTx : GX —)■ it* be the restriction of vr to the orbit GA. 

Definition 62. Define the reduced space = 7rf^{Kfi)/K. 

Remark that <h“^(GA, Kfi)/{G x K) is isomorphic to and that the 

reduced space is non empty if and only if {X, fi) G Gg,k{T*G). 

The [Q, R] = 0 theorem of Meinrenken-Sjamaar relates mG,K{X, fi) to the 
Riemann-Roch number (suitably dehned) of the space As a conse¬ 

quence, we obtain the following theorem. 
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Theorem 63. There exists a decomposition of the cone Cq ^k{T*G) = UaCa, 
in closed solid polyhedral cones Cq such that the following property holds. 

For each a, there exists a non zero quasi-polynomial function pa on Ac © 

Ak such that 

mG,K{\F) = Pa(A,/i) 

*/ (A, £ Ca n (Ag © Ak) ■ 

In particular, for any pair (A, p) of dominant weights contained in Cg,k(T*G), 
the function k mc^K^kX, kp) is of the form: mG,K{kX, kp) = X]i=o 
where Ei{k) are periodic functions of k. This formula is valid for any /c > 0, 
and in particular -E'o(O) = 1. 

An interesting example is the case of K embedded in G = K x K by 
the diagonal. Recall that c\^^ is the multiplicity of 14 in the tensor product 
14 © 14- Thus we obtain 

Corollary 64. Let K embedded in G = K x K by the diagonal. If {X,p, u) 
is in the GG,KiT*G) , the dilated Littlewood-Richarson coefficient k —)■ 
is a quasi-polynomial function of k G {0,1, 2,.. 

To describe the cone Gg,k{T*G) is difficult, and has been the object 
of numerous works, notably Berenstein-Sjamaar, Belkale-Kumar, Kumar, 
Ressayre. We refer to the survey article, [12]. The complete description of 
the multiplicity function mc^K, in particular the decomposition of Gg,k{T*G) 
in UaCa is even more so. However, we will give an algorithm where, given 
as input (A,/i), the output is the dilated function k —)■ mc^K^kX, kp). In 
particular, we can test if the point {X, p) is in the cone Gg,k(T*G) or not, 
according if the output is not zero or zero. 

Consider the set d' C iil of non zero restrictions of the roots A+ to iff. 

We say that is the list of restricted roots. 

Recall that an hyperplane in is ^/-admissible if it is spanned by el¬ 
ements of 4'. The set A = .4,(4^) of admissible hyperplanes is finite. For 
H ^ A, consider JA G h such that H = X^. Let Wg be the Weyl group of 
G. If A G h © tg, consider wX G tg and the hyperplane 

H{w) = {(4 n) G (4 wX) - (z/, W) = 0}. 

We obtain a finite set of hyperplanes IF in R* © ii^. 

Consider a connected component r of the complement of the union of 
the hyperplanes H{w), where H varies over admissible hyperplanes in 
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and w in the Weyl gronp of G, in other words r is a tope for the system of 
hyperplanes J^. So r is an open conic subset of it* © itj. Thus, if (^, z/) e r, 
for any admissible hyperplane H ^ A with equation X, and any w G Wg, we 
have 

{^,wX)-{u,X)^0. ( 7 ) 

The following proposition follows from the description of the Duistermaat- 
Heckman measure |22j . 

Proposition 65. The facets of the cones generates hyperplanes belonging 
to the family T. 

Thus the following lemma follows, and will be useful. 

Lemma 66. Fix a cone Cq. 

• If T is a tope, then r fl it* ^ >o is either contained in Ca, or is disjoint 
from Ca- 

• The closed cone Cq is the union of the closures of the sets r fl it* {>o 
over the r such that r fl Cq is non empty. 

Remark that there might be several topes r needed to obtain Cq. 

We rephrase Theorem [63] as follows. 

Proposition 67 (“Continuity property” ofmc^K)- • If t is a tope, the 
function (A,/i) —)■ is given by a guasi polynomial function 

Ur on r n Ag,k,>0- 

• // r n it*{ >o is contained in the cone Cgxk(T*G), o.nd if (A,p) G 
r n Ag,k,>o, then 

mG,K{\li) = n^(A,/i). 

We hnally state a result on the behavior of the function mG,K{^,T) 
the boundary of the polyhedral cone CGxKiT*G). 

If X G itf, we have an injection K{X) C G{X). Let X E iit, H = X-^, 
and w E W such that {f,wX) — {y,X) > 0 for all [f, v) E Ggxk{T*G). We 
assume that H{w) contains an element (A, /i) G it^ ^ ^q, in other words that 
F is a regular face. Let F = H{w) fl Ggxk{T*G) be the corresponding face 
of the polyhedral cone Ggxk{T*G). 
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As w is defined modulo the Weyl group of the stabilizer of X, we may 
choose w such that, if A is dominant for G, then is dominant for 

G{X). If (A,/i) G F, the couple (wX,^) is a couple of dominant weights for 
{G{X), K{X)). The following proposition follows again from the [Q,R] = 0 
theorem of Meinrenken-Sjamaar (see also USD- 

Proposition 68. For any (A, p) G F, we have mG,K{X, /i) = mG{x),K{x){w\, /i). 

A proof of this theorem using surjectivity of the restriction of homolorphic 
sections, with given invariance, is given in Ressayre m- 

5.2 Branching theorem: a piecewise quasi-polynomial 
formula 

Our aim is to give an explicit quasi-polynomial formula for mG,x(A,/n) on a 
tope r in terms of iterated residues. 

We assume, for simplicity, that iii contains a regular (with respect to A{) 
element X which is regular also for Ag (this is not always the case). We 
use this element to dehne positive compatible root systems A+ and A^ as 
follows: A+ := {a G Ag,Q;(A) > 0} and A^ = {a G A{,a(A) > 0}. Given 
a G A+, denote by a the restriction a\^^. Then 0 7 ^ a G A^. Thus the system 
T in itj consists on the restrictions of A+ repeated with multiplicities: 

T = [a, a G A^]. 

The system T contains A^. By our construction, all elements in T satisfy 
(V^,A)>0. 

We denote by tI/\A^ the list where we have removed A^ from T. More 
precisely, if a G Aj*" occur in T with multiplicity m > 0, then a G Aj*" occur 
in d^\Aj'' with multiplicity m — 1 . 

Example 69. 

Let G = SU{n) and K = SU{ni) x SU{n2) with n = nin 2 . We consider 
t the Cartan subalgebra of g given by the diagonal matrices of trace zero and 
tfc = ti X t 2 the Cartan subalgebras of t given by the corresponding diagonal 
matrices. We take the embedding from iii x R 2 —^ R given by 

diag{ai, ..., a„J x diag{bi,bn^) 
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diag{ai+hi, 02 +^ 1 ,..., a„j+ 6 i, ai+ 62 , • • •, Q,„j+ 62 , • • •, fli+^n 2 ) • • •) ®ni+^n 2 )- 
We take the lexicographic order. The list of restricted roots is thus the list 

= [{ai-Qj+ hk-hij\. ( 8 ) 

There i,j varies between 1 and ni, and k,^ and varies between 1 and 77 , 2 . 
The couple ( 7 ,/c) being different from (j, •^), so all restricted roots are non 
zero. This lexicographic order is compatible, as one check easily, that the 
restrictions are not zero and we get exactly the order Oi > 02 > • • • > on 
ti and similarly for t 2 . 

Explicitly, we can take the embedding of the element 


X = diag{ni ,..., 1) x diag{{n 2 — l)ni + 


For rii = 2,77.2 = 3 then 

X = diag{2, 1), xdiag{5, 3,1) 

and the embedded element is diag{7, 6, 5,4, 3, 2). 

We do not know how to compute the set .4,(T) of admissible hyperplanes 
for the system T, for any ( 77 , 1 , 772 ). We computed it for a few examples (see 
also 0), but we do not see a general pattern. Furthermore, the cardinal of 
the set ^( 47 ) (up to Weyl group action) seems to grow quickly. □ 


Let A G Ag^>o be the highest weight of an irreducible representation of G. 
The character xa of is given by the Hermann Weyl character formula: 


AaItc ^2 


^w(X) 






Restricting on T^, we obtain: 


^w(\) 




llaeA 


by 


Let w G Wg. Let z G (t{)c. Consider the meromorphic function of z given 








55 









Define 



n (1 - 

/3eA+ 




Consider a tope r C it* © it| for the system of hyperplanes Let 
z/) G r. Then for each w G Wg, the element w{^) — z/ is T-regular. Given 
w G Wg, we note a(ta(^) — u) the tope for T, which contains the element 
■*^(0 ~ The tope a{w{^) — z/) depends only from w and r, so we denote it 
by a^. We have dehned the set C>iS(T, a^) of adapted basis to the tope 
(Def. |49]). Let L^^ be the dual lattice to Ak, and let q = g(\l/) be the index 
of T. 


Proposition 70. Let t be a tope in it* © ilj for the system T. Define for 
(A, /i) G Ag © Ak, 

Pr{KL)=Yl Y1 Y1 +— 7 )- 

wjeWj 7Grx/grx o'e05('i',a5,) ^ 

Then, the function Pr is a quasi-polynomial function on Aq © Ak- 


Proof. Explicitly, 


Pr{\,p) = 


ZEE 

-l&'K/qTK a&OS{^,al,) 


Res-^- 


g©(A)-/.,^+^7> 


naGA+(^ 


g-©(a),^+^7>^ 


If (i-« 

/36A+ 



□ 


Theorem 71. Let r be a tope, and let (A,/i) G r H Ag,k,>q- Then 


1. tf{\,fi)^CG,K{T*G), 


mG,K{><,L) = Pr(A,/i). 

2. if {X, fi) G Cg,k{T*G), and the tope r intersect Cg,k{T*G), then 

mG,K{.X,p) = Pr{,X,p). 
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To test if a regular element (A,/i) is in the cone Cg,k(T*G). Take the 
tope r containing (A,p). The point (A,/i) is in the cone if and only if the 
quasi polynomial function k —)■ Pr^kX, kp) is non zero. 

Before going into the proof of this theorem, let us make some remarks of 
how to somewhat reduce the complexity of this formula.. 

Remark 72. 

Since the system T contains the system A^, the function could be 

written as function with the smaller denominator npe 4 '/A+ (1 — e and 

the residues could be taken over adapted OS basis of the system \I/\Aj^. Not 
so many w giving a non zero contribution to the formula. Indeed, at least w 
has to be such that wX — p is in the cone generated by the restricted roots. 
This is the so called valid permutations for (A,/i) dehned by Cochet in [T5] . 
and his algorithm construct them recursively. For example, as computed by 
Pamela Harris, [21], when T is the system f positive roots of A the highest 
root and /i = 0, the number of w such that wX is in the cone generated by 
the positive roots is the Fibonacci number /(r), much smaller that (r + 1)!, 
the order of the Weyl group. 

□ 


Proof. Consider 


Ta|t^ 


E 


n. 


toGWg i 


g«)(A) 

(1 - e-^)’ 


Using our regular element X, we rewrite this formula polarizing the linear 
form w{a): if {w{a),X) < 0, we replace w{a) by its opposite; we then make 
use of the identity 

Precisely, write T = U as the disjoint union of the two sets: 


= {wa, a e Aj", {w{a),X) > 0} 
a e A+, {w{a),X) < 0}. 

If = 1^2 I and obtain that x\\tj^ is equal 

to 


E 


3«'(A) 






E 

UlGWg 
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Lemma 73. The following holds on Tk: 


X\= ^ ^ (-1)""'Pvi>(m;(A)-/ i)e^ 

W'GWg /xeA^ 

where is the partition function determined by the restricted roots So 
’TiG,TK(A,/i) = - p.). (9) 

w&Wf, 

When K is the maximal torus Tg, the formula above is Kostant multi¬ 
plicity formula for a weight [33]. The formula (]9|) is obtained by the same 
method. 

Let us now use the formula 

'mG,K{\p)= 'Yh e(w})mG,Tif(A,/i-M)(pt)+Pt)- 

iZiGWe 

We obtain for (A,/i) G Ag,k,>o 


mG,K(A, h) = ^ (-1 )'*’"Pvi,(m;(A) + g^ - {p - w{pi) + pi)) 

wGWt wGWg 


(we may rewrite this expression as a sum of partition functions for 
obtaining Heckman formula, 1221,H). 

Suppose first that (A, p) G rnAG,_ft:,>o is regular and ’’very” far away from 
all the walls H{w). So for all w, w, the element w{\) — p + gw + — pt 

is still in the tope for T containing w{\) — p. We thus can employ the 
iterated residue formula (Formula[6]) for V^{w{X) + Pw — ~ 'Mi(Pf) + Pt)) on 

We obtain that mG,K{X,p) is equal to 

{w{\)+g^-{ti-w{pt)+pt),z+if-'y) 

y e(h))(-i)'’“ y y y Res^- --—5——. 

W&^l «)GW0 7er-f^/gr^ fTGG<S('I',0j,) llpG'I'V-'- ® ) 


Inverting the polarization process, we rewrite 


(-ir 


g(«)(A)+g^,z+^7> 


g(«;(A),z+^7> 






n 


osa; 


(1 -e 


— {wa,z+^^'^) 
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So, remembering that n^eA+ we obtain 

mG,K{KlJ^) = Pr(A,/U) 

when (A,/i) G r fl Ag^k,>o is regular and ’’very” far away from all the walls 
H{w). Theorem [HS] asserts that is given by a quasi polynomial 

formula on r Pi Ag,k,>o or even on r fl Ag,k,>o if t fl ft* ^ >o is contained in 
Cg,k{T*G). 

Using Lemma [28] we obtain Theorem [TO] □ 

5.3 Singular case 

When the stabilizer of A is large, that is (A, Hq) is equal to 0 for a large 
number of a, then we can rewrite the Hermann Weyl formula for the character 
in a way that takes advantage of this. 

Fix a subset S of the simple roots of A+, and let 1 be the Levi subalgebra 
of 0 , with simple root system E. Let be its positive root system. 

Let ft* the set of the elements ^ G ft* such that {^, Ha) = 0 for all 
a G E. We dehne consistently Aq = Ag n ft* s, a lattice in ft* 

^G,K,>Q = ^G,K,>0 n ft* 2 , and similarly. 

We also dehne 

CE,KiT*G) = {(e,^) e Gg,k{T*G)-^ G 
In other words, u must belong to the projection on ft* of the singular orbit 

Ge 

The cone Gqj^{T*G) is solid in ft* s, if and only if there exists ^ G ft* ^ 
such that the projection on ft* of the singular orbit G^ has a non zero interior 
in ft*. In other words the Kirwan polytope n{G^) fl ftj is solid. 

Example 74. 

Consider the embedding of iC = U{n 2 ) x U{n^)/Z in G = U{n 2 n^)- Here 
A 2 and Z is the subgroup {t 2 ld,t^Id} of the center of U{n 2 ) x U{n^) 
with ^ 2^3 = 1- We take A a weight of G with more than two non zero 
coordinates. Then 7 r(GA) has interior in ft*. 

□ 

Let 

Au = A+\A+ 
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and denote by >V( the Weyl gronp of [. 

For any A G >o H 5 ., we can write the character formula on Tc and 
the restriction on Tk as: 


'X.Atq 




^w{X) 


«)GWg/W| 






E 

tOGVVg/Wl 


^w(\) 


In the regular case, the Levi component is just the Cartan subalgebra , the 
parabolic is the Borel subalgebra, and Au = A+. 

To compute mG,K{^, jj) for AG,>onitj by iterated residues, it is tempting 
to replace the function 5'“^ by the function 




Yl,w , 
X,fi ' 


z = 


n(i 






/ 36 A+ 




We consider the system of hyperplanes in it* © it^ defined by the 
equations (^, wX) — (n, X) = 0, where X is an equation for a d'-admissible 
hyperplane, and w G Wg. If (^, u) G it* ^ © itj is in a tope Ts for then 
(^, u) is in a unique tope r for and determines a tope ajj in itj. 

The following proposition is clear. 

Proposition 75. Let 

p^(A,p)= ^es^‘^A,r(^+^)- 

W£WJV\>1 'idVK/qVK (T&OS{^,al,) ^ 


Then p^(A, p) is a quasi polynomial function on Aq © Ak- 

Remark that if ts is a tope for X-£, then rs H Cqj^{T*G) is empty if 
Cq is solid. If the cone Cq i® solid, it is the union of the 

closures of the open sets fl t^contained in ;^(T*G). 

Theorem 76. Let ts a tope init* 5 ^©it| forXj], and let (A,p) G TsPiA^^^ ^q. 
Then 
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1 . 


mG,x(A,/i) =p^(A,/i). 

2. if {\, fi) G Cq j^(T*G), and the tope intersect Gq j^(T*G), then 

mG,x(A,/i) = p^(A,/i). 


Proof. If the set n Gq j^{T*G), is contained in r fl Gg,k{T*G), so we 
know (by before) that on r fl the function mG,K is given by a quasi 

polynomial formula, so a fortiori its restriction to ts H ^ >q. So it is 
sufficient to prove that when fl A^ >q is sufficiently far away from all 
walls belonging to then mG,K{^,p) coincide with p^. 

The proof is very similar to the preceding proof, so we skip details. 

Let w G Wg. Because of our assumption on compatible systems, namely 
the existence of a regular element X, we can dehne ^ U ^ with 

= {wa,a G Au, («;(«),X) > 0}, = {-w{a),a G A^., {w{a),X) < 

0}. Elements in \ku,,u are positive on X, so is contained in T. In contrast 
to the regular case, \ku,,u depends of w and may not contain A^. 

Let J and {TR,x)<o then we obtain 


AA|t^ 


E 

WGWgfWi 




and 


mG,K(\p) = e(w) (-iy^'P^^,u(wW+gE- (M-w(pt) + pt)). 

weWt WGWg/Wl 


The point (A, p) being in te, the point to(A) — /r is in We can assume 
that (A,/i) is sufficiently far away from all walls, so that w(A) + ~ (P ~ 

w(pt) + pt)) is also in aj),, so that we can apply the iterated residue formula 
for 

L^'i'^,u(w^(A) + g^-(M- w(pi) + pe)) 

as a sum of iterated residue with respect to adapted OS basis for Then 
we proceed as in the preceding case, reversing the polarization process, and 
obtain 

mG,K(A,p) =pf(A,p) 
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provided (A,p) is in and snfficiently for away from the walls. 

If r n Cq j^{T*G) C r n Cg,k{T*G) the formula is quasi polynomial, so 
we obtain our theorem. 

□ 


6 Kronecker coefficients and examples 

We describe now our approach to compute Kronecker coefficients, the result 
is summarized in Corollary 1771 

Consider N = rii- ■ - Ug and assume that ni is the maximum of the n*. 
Write T-i = and consider the action of U (ni) x • • • x f/(n^) 

in Sym{'H). 

Thus 


Symin) = Y ,® ® ® ® 

We want to compute g{yi ,..., ys) and the dilated coefficients. 

Let M = n 2 ■ ■ - Ug, and consider the embedding oi K = U ( 77 , 2 ) x • • • xU{ng) 
in G = U (M) determined by: (^ 2 , • • • , kg){v 2 ■ ,Vg) = ^ 2^2 <8) • • • kgVg, as 
we explained in Example [5l 

Using the Cauchy formula [HI we can write 

%m(C^) = 0 C^) = Y 

>0 


Write the decomposition of 


restricted to K: 


yU{M) 


V2G[/(n2),- ,/ise;7K)"^G,A:(hn h2 ® 




yU{ns) 

fj,s 


Remember that the polynomial irreducible representation of U{nk) are pa¬ 
rameterized by the highest weight 7 = [ 71 ,... , 7 ^*,] with 71 > 72 > • • • > 
Irik ^ 0 and that I 7 I = Yhili- Taking care of the fact that if pi 0 • • • 0 /i^ 
occur in Sym{C^^ 0 ■ ■ ■ <E> C”"*) we must have that the actions on the centers 
must be the same, we obtain: 


Corollary 77. 


g(iui 0---0yg) 


mG,K{yi,y20 ■ ■ ■ 0 g^s) if \gi\ = \g2 
0 otherwise 


Ihs 
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We take advantage of the formula to somewhat reduce the computation 
of a tensor product with s factors to an analogous computation with s — 1 
factors. So we will compute some of the Kronecker coefficients as a corollary 
of the branching theorem. 

The case 0 (8) corresponds to give explicit formulae for the 

decomposition of the representation of the symmetric group associated to 
partitions with at most two rows. Complete expressions for these functions 
have already been obtained by Briand et al. 

We list the examples we can compute in Subsection 16.21 For example, we 
can compute the dilated Kronecker coefficients for 0 0 C^, as well as 

some examples for 0 0 C"^. 

6.1 The algorithm to compute Kronecker coefficients 

We refer to Section 15.31 and Section |6] for the notation. 

We are given a sequence of s strictly positive integers [ui,..., Ug] and for 
each integer rii a sequence z/j of integers: Vi = - ■ ■ , r'^,] with i'\ > i' 2 > 

• • • > > 0. Each Vi parametrizes an irreducible polynomial representation 

of Uijii) of highest weight Vi. Write N = ni=i M = ni =2 want 

to compute the Kronecker coefficients g{kvi, • • • , kvg) dilated by an integer k 
that is the multiplicity of the tensor product representation kvi 0 • • • 0 kvg in 
Sym{C^). Our approach uses Cauchy formula together with the computation 
of the branching coefficients to reduce the number of parameters. We may 
assume that rii < M and that \vi\ = \v 2 \ = ■ ■ ■ = \vs\. 

We set G = U{M) and K = U{n 2 ) x • • • x f/(n,). 

The first reduction step is: 

• If \vi \ = \v 2 \ = ■ ■■ = \i^s\ then g{vi 0 • • - 0 z/^) = mG,x(d'i, V 2 ®---®^s) 
where vi is the highest weight representation of U (M) obtained by v adding 
M — rii zeros and the branching coefficient mG,K is computed in Theorem [76] 
via the function defined in Proposition [TS] 

Let us write A = z7]^, /i = z/ 20 - • •0z/s. If ni < M, A is a singular weight for 
the group U (M). Denote by S the set of simple roots [e „^+2 — ^m+i, ■ ■ ■ — 

cm-i] of U{M). Let I = u(M — rii) the Lie algebra with this simple root 
system. We have (A, Ha) = 0 for all a G S. 

Then (A, p) G Aq >q © Ak>o with the notations as in 15.31 Let us review 
the key steps of the algorithm to compute mc^K- See the discussion in 14.9.81 
outlining the limits of the implementation. 
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Given as input (A,/i), we wish to compute the branching coefficients. 
Recall that : 

mG,K{Kf^)= +( 10 ) 

w€Wg/yVt'yerK/qrK<TGOS{^,ai,) ^ 

One of the tricky point in computing the right hand side of equation flTOl) is 
to find a J^s-fope rs such that (A, /i) G T^. We do this by computing a regular 
point inside the Kirwan cone sufficiently small and deform (A,/i) along the 
line from (A,/i) to this interior point. 

1. We list all the equations X of the X- admissible hyperplanes. 

2. For each such equation given by X and for w G Wg, we compute 
H{X,w) = {X,wX) — {fi,X). As (A,/i) is in the lattice of weights, 
and X in the dual lattice, H{X,w) is an integer. 

Remember A,p are our hxed input. 

3. If H (X, w) 0 \/w, X, then (A, /i) is Xs-regular and then it is in a tope 
Ts. A fortiori it is in a unique X tope r and therefore w{X) — /i is in a 
unique tope G 

In conclusion we can compute a^). 

4. Else if H{w,X) = 0 for some X and w, then we deform as follows: 

(a) We hnd e = (ei,e 2 ) in the interior of . We can End this 

point, in the cases we treat, because either we have the equations 
of the Kirwan cone, either we know some points in the Kirwan 
cone by directly computing projections. 

(b) We rescale e so that |(taX, ei) — (X, € 2 )! < 1 / 2 , so (A,/i) +te stays 
in the same tope for all 0 < t < 1 . 

(c) We define (A^e/, l^-def) as (A + Ci, + 62 ). 

5. We can now pick the tope r defined by (A^e/, A^de/) and compute OS (\k, ajj,) 
as in step |3l 
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6 . Now we compute 


sf::{z) = 


n(i- 

/3eA+ 


,-(/ 3,2 




^{w{\)-ll,z) 




and the residue along an OS basis with an appropriate series expansion. 

7. Finally to compute mc,K{^, we have to sum the contribution from 
w G Wg/W(, over the set 7 G Tk/q^k and a G OS{'^,al^). Each 
individual term of these two sums, that is if we £x a 7 and a a, is 
easy to compute, but there can be really many of these terms, making 
possibly the computation very long. 


Observe once again that in particular, we can test if the point (A, p) is in the 
cone Cq j^{T*G) or not, according if the output is not zero or zero. 

It is not more difhcult to compute with symbolic variables (A,/i) be¬ 
longing to the closure of a tope. However to describe all possible topes 
(the chamber decomposition of Cx('H)) seems very difficult. So our input is 
(A,/i), the output is either the numeric value, either the dilated coefficient 
k —)■ mc^K^kX, k^), or (in low dimensions), a tope r containing A,/i in its 
closure and the quasi-polynomial function in both variables A,p coinciding 
with mG,K{X,fi) on the tope r. 


Remark 78. Rectangular tableaux 


Remark that if pi is a rectangular tableau, then A is even more singular. 
This enable us to compute more easily using a larger set S (reducing then 
the number of parabolic roots). When all /ij are rectangular tableaux, this 
corresponds to the case of Hilbert series. We list the corresponding results 
in the last subsection 16.31 


6.2 Examples 

The hrst two examples have already been treated in the literature. 
Example 79. The case of 3-qbits : ( 8 ) 0 

This case has been treated in complete details in 14.61 and is due to [TT] . 
Example 80. The case of : 0 
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This example has been studied in complete details by m- The number 
of chambers of polynomiality is 74 and on each chamber the quasipolynomial 
is of degree 2 . 

The multiplicity function k —)■ m{kX, /c/i, ku) is a quasi polynomial func¬ 
tion of the form 

/(i) + (-l)‘c 

where f{k) is a polynomial of at most degree two and c is a constant. Here 
is an example. For A = [5, 3, 2,1], p = z/ = [ 6 , 5], then we obtain 

1/4 P + l/2k + 5/8 + 3/8 (-1)^ 


Remark that all points (a,/?, 7 ) in the boundary of the Kirwan cone are 
stable, thus g{ka, kf3, k'j) is 0 or 1. 

Example 81. The case of 4-Qbits : (8) <8 <8) 


We consider the action of of iF = U{2) x 7/(2) x 7/(2) x 7/(2) on 77 = 

00 ( 8 ) with (fci,/c 2 , ^ 3 , ^ 4 ) acting as/i;i 0 /c 2 0/^3 0 ^ 4 . The Kirwan 
polytope has been described by Higuchi-Sudbery-Szulc, [23] . 

We have no idea of the number chambers for polynomiality. Neverthe¬ 
less, given highest weights a,{3,'y,6, we can compute g{ka, k(3, k'j, k6) as a 
periodic polynomial in k. 

It is a polynomial of degree at most 7 and period 6 . Here is an example. 
When a = /3 = j = 6= [2,1] then we compute: 


m{ka, kjS, kj, k6) = 


23 


,7 13 ,e 155 ,5 139 ,4 

k^ + ——k^ H- 


241920 5760 6912 

9799 , 5 

-h (-1) - 

11520 ^ ^ 256 


r + 


1152 
738545 
V 32256 


7 81601 1 , 

--^-(- 1 )^ 

V207360 1536^ ^ J 

179 A 

+ (-1)^- k + P{k) 

^ ’ 1536 / ^ ’ 


where 


P{k) = 





1 A 5279 51 

- 9\d ^ -^- 

243 J 6912 256 


(- 1 )' 


is of period 6 and 6 * is a primitive root 6^ = 1. 
0,1, 2, 3,4, 5 are 

5725 76 77 77 5597 
10368’ 10368 


The values of P{k) on 
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Figure [7] shows the Duistermaat-Heckman measure for four qbits. The 
drawing is along the line from = [[1/2,1/2], [1/2,1/2], [1/2,1/2], [1/2,1/2]] 
to Vtop = [[1,0], [1,0], [1,0], [1,0]]. The function t DH'^{tVmin +-t)vtop) 
is a spline of degree 7 with singularities at f = [0,1,1,1,1]. 



Figure 7: Duistermaat-Heckman measure for four qbits 

Example 82. The case of : C® <8 <8) 

When n 2 = 3, = 2, it is sufficient to consider the case when rii = 6. 

Then the multiplicity function k —)■ m{kX, /c/r, ku) is a quasi polynomial 
function of the form 

f{k) + (-l)'‘g(k) + h(k) 

where f{k) is a polynomial of k generally of degree 8, g{k) of degree 2 and 
h{k) is a periodic function of A; mod 6. 

Here is an example where the degree of the polynome is the expected one. 
We £x A = [1500,1052,940,492,268,156], ja = [2110,1438,860], i/ = 
[2748,1660] and compute: 


20160143036868818273 


540 

311654396584249 


k^ 


m{k\, kfi, ki/) = 
2401100429169038668 


270 




945 

11934644414969 




236265968398572733 


3240 




232 


1080 

— + — e 

81 243 


,, 53468894201 ,, 44636639 , „ 
k^-\ -Ph--F 




810 

1243 


180 


378739 

630 


k+ 


243 
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where 0 is a third primitive root of 1. The term of the degree zero in 
m{kX, kjj, kv) is a periodic function of k whose values are 

25 76 ^ 25 76' 

We hnish by noticing that for fc = 1... 6 , the values of m{kX, kfi, ku) are 
given by the following list 

[1,39948532219001323,9887333657571493818,250556011548476811713, 

2488623801870416780185,14783083490287618355455], 

We now compute examples of multiplicity on the walls. 

The walls of the Kirwan cone have been described by Klyachko in [32] , 
Given A = [Ai, A 2 , A 3 , A 4 , A 5 , Ag], p = [/ii,/i 2 ,/is] and u = [iii,z/ 2 ,], the walls 
are by the following 5 types of inequalities in A, fi, v. 

More precisely, for each of the inequality below, there is a particular 
non empty subset S' of Eg x S 3 x S 2 (where S^ is the of permutations of k 
elements) computed by Klyachko such that the permuted inequality is a wall 
of the corresponding Kirwan cone. 


type 

equation 

I 

i/i — 1/2 ~ Ai — A 2 — A 3 + A 4 + As + Ag < 0 

II 

/Ti + /i 2 ~ 2 /r 3 — Ai — A 2 — A 3 — A 4 + 2 As + 2 Ag < 0 

III 

2 /ii — 2 /r 3 + z/i — z /2 ~ 3 Ai — A 2 — A 3 + A 4 + As + 3 Ag < 0 

IV 

2 /ii + 2 /i 2 — 4 /i 3 + 3 1^1 — 3 z /2 — 5 Ai — 5 A 2 + A 3 + A 4 + As + 7 Ag < 0 

V 

4 /ii — 2 /i 2 — 2 /i 3 + 3 1^1 — 3 z /2 — 7 Ai — A 2 — A 3 — A 4 + 5 As + 5 Ag < 0 


The following table list elements v = [A, /i, v] in the relative interior of the 
corresponding wall type facet, together with the value of m{kX, kfi, ku). As 
asserted in Lemma ITTl for the listed value m{kX, kji, ku) has the maximum 
possible degree. 


vector 

[A,//,i/] 

m{kX, k^, kv) 

VI 

[[288,192,174,120,30, 6], [343, 270,197], [654,156]] 

1 + 17A: 

VII 

[[600,372, 300,156,96,12], [876,552,108], [930, 606]] 

1 + h + + 242154 

VIII 

[[188,140, 92, 52, 20,4], [304,152,40], [340,156]] 

1 

viv 

[[276, 204,120,66, 30,6], [351, 273, 78], [552,150]] 

1 + 36A: 

Vy 

[[276,198,126,66,48,6], [406, 201,113], [536,184]] 

1 + 41A: 
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In particular, we see that all facets of type vjjj consists of stable elements. 

There are 20 such facets (obtained by considering the special permutations com¬ 
puted by Klyachko). 

As we noticed just after Theoreml.Sblwe don’t know how to compute the degree 
when /r is singular, see Subsection 16.31 for the case of three rectangular tableaux. 

Here is an example for which the degree is smaller. Consider A = [9, 7, 5,3, 2,1], 

/i = [9,9,9], u = [14,13], then 

m{kX, kn, ku) is given by the following formula 

(i (-1)'= + i) k+% k^+W2 (-2 G + 8)+^ (2 0 + 10) + 

JA 4 - A 

144 16 1 1 

Here 6 is again a primitive root 6^ = 1. Thus the term of degree zero is a 
periodic function r oi k such that 

[r(0), r(l), r(2), r(3), r(4), r(5)] = [1, 71/216,17/27,5/8,19/27, 55/216] 

Of course, the value of m(0,0,0) is equal to 1. Here m(A, p, i^) = 5 and for instance 
the value m(17A, 17/r, 17z^) = 344715. 

Example 83. The case of 3-quthrits. : O O 

Here TL = C^(8)C^(8)C^. The multiplicity function k m{k\, kpL, kv) is a quasi 
polynomial function of the form 

f{k) + {ifg{k) + h{k) 

where f{k) is a polynomial of k of degree at most 11 and h{k) is a periodic function 
of k (mod 12). The actual numerical values are computed in a rather quick time. 

Let us give an example of the dilated Kronecker coefficient. The periodic term 
for the coefficient of degree 0 m{kX, kg,, ku) with X = g = v = [A,2>,2] is given by 

1166651 13403 29899 59 1166651 235 980027 59 32203 13403 980027 ' 
5308416’ 20736’ 65536’^’ 5308416’^’ 5308416’^’ 65536’ 20736’ 5308416 

In this case m{kX,kg,kv) has precisely degree 11. 

Here is a numerical example. For A = [5, 2,1], = [4,2,2], v = [3,3,2] we 

compute m(A, g, v) = 4. 

We close with one last example that we had promised. 

Example 84. 
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Here V. = 0 C^. 


9([M1. 1,1], Ml. 1.1], t[2, 111) = j (-1)‘ + ^ + (1 - 9) + (2 + 9) + ifc 

and for instance the list of the values for 0 < A: < 30 are 


[1,0,1,1,1,1, 2,1,2, 2, 2, 2, 3, 2, 3, 3, 3, 3,4,3,4,4,4, 4, 5, 4, 5, 5, 5, 5, 6] 


So the saturation factor is 2. 


6.3 Rectangular tableaux and Hilbert series 

We give a list of the Kronecker coefficients for the following situation of rectangular 
tableaux. We use the following notations: (C^)^ = (g) ® C^, [[1,1]]^ = 

[[1,1], [1,1], [1,1]], (8) and similarly. In the following Table 

the second column refers to the choice of the parameters [A, p, ly] and the third 
column to the value of the Hilbert series m{k)t^ where m{k) is the Kronecher 
coefficient g{kX,kfi,ki'). 


type 

parameters 

value 

(C2)® 

[|i.iii“ 

1 

(C2)4 

CT 

I 

(C2)^ 

CT 

HS 22222 

(C3)® 

[[MMl? 


C[4>3>3] 

[[3,3,3,3], [4,4,4], [4,4,4]] 

1+P 


where 

HS 22222 = 1], 1], ^[1, l],k[l, l]])t'= = P{t)/Q{t) 


where 


P{t) = + 16+ 9+ 82+ 145 + 383+ 770 1‘^^+ 

1659+3024+5604+9664 P +15594 P+24659 P+36611+52409+ 
71847+95014+119947+146849+ 172742+195358+214238+ 
2256991^^+2297521^^+2256991^^+2142381^^+1953581^3+17274212^+1468491^^ + 

119947 i2°+95014 1^^+718471^3 _^524091^^+366111^®+246591^^+155941^'‘+96641^3 

56041^2 30241^1 + 16591^° + 7701^ + 3831® + 1451'^ + 821® + 91® + 161^ + 1 
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and 


Q{t) = (1 - t2)5(i _ _ ^5)^^ _ ^6)5^ 

We remark that for the 5-qubits the result in [38] correspond to the series 
and has a misprint on the value of the coefficient for n = 42 
(corresponding to the coefficient of in our formula), as the numerator is not 
palindromic. So the value o„ for n = 42 in [38] has to be replaced by 146849. 

We report for completeness the value of the Kronecker coefficients in the exam¬ 
ples considered, we omit the actual expression for the Kronecker coefficients in the 
5-qbits case and the one for 5(A:[3,3, 3,3], A:[4,4,4], A;[4,4,4]) because the formula 
is too long to be reproduced here. 




<^o 11 1 

(2 + e) + —02^ (1 - 0) + 

36 4 ^ ^ 27 ^ 27 ^ ' 


29 1 

48 16 


+ — (-1) U + —+ 


16 
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5(A:[1,1,1],A:[1,1,1],A:[1,1,1]) = 


107 9 


288 32 


^ (- 1 )'" + (1 + (- 1 )^ ) + (1 + (- 1 )''^') 


\ A:+l 


16 




’k+1 


16 


9 


9 


16 


16 


48 


where 0 is a third root of unity. For g{k\l, 1, 1], k[l, 1, 1], A;[l, 1,1]) we report, as an 
example, the expressions on cosets. We have twelve cosets and thus a sequence of 
12 polynomials given by the following list 

r . 1 , 1,2 7 1 , 1,2 5 1 , 1,2 7 1 , 1,2 

^ 4 48 48 8 48 12 4 48 16 8 48 

2 1 , 1,2 7 1 , 1,2 3 1 , 1,2 5 1 , 1,2 

3 4 48 48 8 48 4 4 48 48 8 48 

2 1, 1,03 1, l,o5 1, l,o5 1, 1,0, 

3^ 4 48 ’ 16 8 48 ’ 12 4 48 ’ 48 8 48 ^ 


The following is the list of values of the Kronecker coefficients computed by 
the above formula for 0 < A: < 20 : 


[1,0,1,1, 2,1, 3,2,4, 3, 5,4, 7, 5, 8, 7,10,8,12,10,14]. 


Once again the saturation factor is 2. 
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